AD-A058  495 

UNCLASSIFIED 


OHIO  STATE  UNI V COLUMBUS  ELECTROSCIENCE  LAB  F/G  9/5 

A HYBRID  TECHNIQUE  FOR  COMBINING  THE  MOMENT  METHOD  TREATMENT  OF — ETC(U) 
JUL  78  E P EKELMANf  G A THIELE  N00014-76-C-0573 

ESL-784372-5  NL 


1 of  rtf 
AD 

A058495 

■ 

1 I 

J 

a 

Pi 



J 

COMBINING  THE  MOMENT  METHOD  TREATMENT 
WITH  TOE  GTD  FOR  CURVED  SURFACES 

Ernest  P.  Ekelman,  Jr. 

Gary  A.  Thiele 


The  Ohio  Stale  University 

ElectroScimce  Laboratory 

. 4 ITf  A .1  4 irt  m fm  1 

IMpariinvill  or  ctecincoi  engineering 

Columbus,  Ohio  43212 


TECHNICAL  REPORT  784372-5 
Contract  N00014-76-C-0573 
July  1978 


JSfe 


k-  4. "'-s- ’’ sri  5 


mmm  |l 

1 t/ff-vJ 

: . >V 


<••  i-  :.  .rl  ,■/  f ;y.  ./  <•  /■  :V  .•  • uSt&Zjj!  :? ' <, 

E if IwPp 


Department  of  the  Navy 
Office  of  Naval  Research 
800  N.  Quincy  Street 
Arlington,  Virginia  22217 


$*!$?■ 


78  08 


NOTICES 


Wien  Government  drawings*  specifications*  or  other  data  are 
used  for  any  purpose  other  than  In  connection  with  a definitely 
related  Government  procurement  operation*  the  United  States 
Government  thereby  Incurs  no  responsibility  nor  any  obligation 
whatsoever,  and  the  fact  that  the  Government  may  have  formulated, 
furnished,  or  in  any  way  supplied  the  said  drawings,  specifications, 
or  other  data*  Is  not  to  be  regarded  by  Implication  or  otherwise  as 
In  any  manner  licensing  the  bolder  or  any  other  person  or  corporation 
«r  conveying  any  rights  or  permission  to  manufacture,  use,  or  sell 
any  patented  Invention  that  May  In  any  way  be  related  thereto. 


UNCLASSIFIED 

SECURITY  CLASSIFICATION  OF  THIS  PACE  (Whan  Data  Bntarad) 


REPORT  DOCUMENTATION  PAGE 


READ  INSTRUCTIONS 
BEFORE  COMPLETING  FORM 


1.  REPORT  NUMBER 


2.  GOVT  ACCESSION  NO 


3.  RECIPIENT'S  CATALOG  NUMBER 


TITLE  (and  Subtllla) 


HYBRID  TECHNIQUE  FOR  CQMBINING  THE  JJOMENT 
THOD  TREATMENT  OF  WIRTANTENNAS  WITIT THE  / 
D FOR  gjRVED  SURFACES,  ^ ' 


V*.  TYPF  nr  pppnflT  A.amnn  riwFPFn 

Mfjechnical  Report / 


At^HORf^ 


Ernest  P.|Ekelman,  Jr. 
Gary  A. 1 Thiele 


9 


ERFOBMIHlI  ORGAMI7  ATIOM  N AMS-ANO(  ADDRESS  / 

The  Ohio  State  University  ElectroScience 
Laboratory,  Department  of  Electrical 
Engineering,  Columbus,  Ohio  43212 


6 rirnrsntmiiii  nnn  rtrrnpT 

SL— 784372-5  J / 


NUMBER 


CONTRJ 


3R  grant  numberi<; 


Contrki^NOOQl  4-76-C-0573  j 


10.  PROGRAM  ELEMENT.  PROJECT.  TASK 
AREA  6 WORK  UNIT  NUMBERS 


Project  No.  NR  371-108/12- 
29-77  (427) 


ii. 


CONTROLLING  OFFICE  NAME  AND  AOORESS 

Department  of  the  Navy,  Office  of  Naval  Resea 
800  N.  Quincy  Street 
Arlington,  Virginia  22217 


112.  REPORT  DATE 

[July*978  S' 


>3.  NUMBER  OF  PAGES 


126 


I*.  MONITORING  AGENCY  NAME  ft  AODRESSfff  dlttarant  from  Controlling  Otllca) 


IS.  SECURITY  CLASS,  (ol  thla  raporl) 

Unclassified 


is  a.  DECL  ASSIFIC  ATION/ DOWNGRADING 
SCHEDULE 


16.  DISTRIBUTION  STATEMENT  (ol  thla  Raporl) 


This  document  has  been ' approved" 
for  public  release  and  solo;  iia 
distribution  is  unlimited. 


U)_L^p. 


17.  DISTRIBUTION  STATEMENT  ( of  the  abatract  entered  In  Block  20,  it  dlttarant  from  Report) 


18.  SUPPLEMENTARY  NOTES 


The  material  contained  in  this  report  is  also  used  as  a dissertation 
submitted  to  the  Department  of  Electrical  Engineering,  The  Ohio  State 
University  as  partial  fulfillment  for  the  degree  Doctor  of  Philosophy. 


19.  KEY  WORDS  (Continue  on  ravaraa  aida  //  nacaaamry  and  Identity  by  block  number) 


Hybrid  method 
Moment  method 
GTD 


Scattering 

Antennas 

Geometrical  optics 
Diffraction 


ABSTRACT  (Continua  on  ravaraa  aide  It  nacaaamry  and  Identity  by  block  number) 


^ The  hybrid  technique  presented  in  this  paper  is  a method  for 
solving  electromagnetic  problems  in  which  an  antenna  or  other  dis- 
continuity is  located  on  or  near  a conducting  body,  such  as  antennas 
on  ships  or  aircraft.  The  technique  solves  these  kinds  of  problems 
by  properly  analyzing  the  interaction  between  the  antenna  or  scatterer 
and  the  conducting  body.  The  hybrid  technique  accomplishes  this  by 
casting  the  antenna  structure  in  a moment  method  format  then  modifying. 
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20. 

^that  format  to  account  for  the  effects  of  the  conducting  body  via  the 
geometrical  theory  of  diffraction.  The  technique  extends  the  moment 
method  to  handle  many  problems  that  cannot  be  solved  by  GTD  or  the 
moment  method  alone. 

In  general,  arbitrary  radiators  located  on  or  near  canonical 
shapes  or  combinations  thereof  can  be  solved  using  the  hybrid  technique. 
Electromagnetic  parameters  for  which  the  hybrid  technique  can  solve 
include  the  near  and  far  fields,  current  distributions,  impedances, 
and  scattering  data.  In  this  paper,  wire  antennas  are  analyzed  to  find 
their  input  impedance  when  they  are  located  near  perfectly-conducting 
circular  cylinders. 

The  purpose  of  this  paper  is  to  present  the  technique  and 
demonstrate  some  of  its  facility  and  its  accuracy.^.  Three  orthogonal 
orientations  were  identified  and  antennas  to  matctvthem  were  used. 

For  each  case,  the  hybrid  solution  is  checked  with  one  of  three 
independent  solutions;  an  ^-eigenfunction  solution,  image  theory, 
or  experimental  measurement.  The  hybrid  technique  was\  verified  for 
these  orthogonal  orientations.  Thus  a more  general  radiator  with  an 
arbitrary  combination  of  orientations  could  be  solved  v^ith  confidence. 
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In  recent  years  the  analysis  of  electromagnetic  problems  has 
been  advanced  tremendously  by  the  digital  computer  and  use  of  two 
powerful  theories,  the  method  of  moments  and  the  geometrical  theory 
of  diffraction.  The  moment  method,  MM,  is  a numerical  technique  which 
converts  an  integral  equation,  containing  field  terms  and  appropriate 
source  terms  that  describe  the  electromagnetic  problem,  into  a system 
of  linear  algebraic  equations.  The  equations  have  unknown  coefficients 
related  to  an  expansion  of  a current  distribution  on  the  structure  of 
interest.  The  moment  method  is  characterized  as  a low  frequency 
method  since  its  practical  use  is  limited  to  structures  that  are  not 
electrically  large.  The  geometrical  theory  of  diffraction,  GTD,  is  a 
ray  optical  method  whose  solution  is  based  on  the  asymptotic  approxi- 
mation to  the  integral  defining  the  electromagnetic  problem.  GTD  is 
characterized  as  a high  frequency  method  as  it  is  applicable  to 
bodies  that  are  arbitrarily  large  electrically.  Both  of  these  methods 
are  powerful  computational  tools  which  permit  application  to  a wide 
range  of  problems.  In  this  paper  these  two  techniques  are  combined 
into  a single  technique  called  the  "hybrid  technique."  The  hybrid 
technique  permits  one  to  solve  many  problems  that  are  not  solvable  by 
either  technique  alone. 

The  hybrid  technique  presented  in  this  paper  is  a method  for 
solving  electromagnetic  problems  in  which  an  antenna  or  other 
discontinuity  is  located  on  or  near  a conducting  body,  such  as  antennas 
on  ships  or  aircraft,  feed  antennas  near  the  reflecting  surfaces  of 
reflector  antennas,  and  slots  or  other  discontinuities  on  conducting 
surfaces.  The  technique  solves  these  kinds  of  problems  by  properly 
analyzing  the  interaction  between  the  antenna  or  scatterer  and  the 
conducting  body.  The  hybrid  technique  accomplishes  this  by  casting 
the  antenna  structure  in  a moment  method  format  then  modifying  that 
format  to  account  for  the  effects  of  the  conducting  body  via  the 
geometrical  theory  of  diffraction.  The  technique  extends  the 
moment  method  to  handle  many  problems  that  cannot  be  solved  by  GTD 
or  the  moment  method  alone. 

The  basic  hybrid  technique  used  in  this  paper  was  first  described 
in  the  literature  by  Thiele  and  Newhouse  [1].  In  that  paper  the 
technique  was  applied  to  antennas  on  and  near  finite  planar  surfaces. 
The  moment  method  solution  was  modified  to  account  for  the  finite 
planar  surfaces  using  wedge  diffraction  theory.  In  the  present  paper 
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a hybrid  technique  for  combining  the  moment  method  treatment  of  wire 
antennas  with  the  GTD  for  curved  surfaces  will  be  presented.  Specifi- 
cally, wire  antennas  will  be  analyzed  near  perfectly  conducting  circular 
cyl inders. 


The  wire  antenna  will  be  cast  in  a moment  method  format  given  by 
Richmond  [2,3]  in  his  computer  program  for  thin-wire  structures.  The 
use  of  piecewise  sinusoidal  basis  functions  in  this  method  will  be 
exploited.  The  effects  of  the  cylinder  will  be  found  using  GTD.  The 
field  expressions  are  obtained  using  the  uniform  geometrical  theory  of 
diffraction  for  an  edge  given  by  Kouyoumjian  and  Pathak  [4],  The 
application  of  the  GTD  to  the  cylinder  is  based  on  techniques  presented 
by  Marhefka  [5]. 


The  hybrid  technique  as  presented  here  is  applicable  to  a much 
broader  class  of  problems  than  those  that  are  demonstrated.  In  general, 
arbitrary  radiators  located  on  or  near  canonical  shapes  or  combinations 
thereof  can  be  solved  using  the  hybrid  technique.  The  arbitrary 
radiator  would  be  modeled  with  a thin-wire  grid  and  set  up  in  the 
moment  method  format.  The  conducting  body  is  restricted  to  canonical 
shapes  for  which  a GTD  solution  exists.  Electromagnetic  parameters 
for  which  the  hybrid  technique  can  solve  include  the  near  and  far 
fields,  current  distributions,  impedances,  and  scattering  data.  Of 
course,  many  geometries  present  special  problems  for  which  the 
technique  will  fail  and  therefore  need  to  be  modified.  For  example, 
part  of  the  radiator  might  lie  in  the  caustic  of  a needed  GTD  field. 


It  is  impossible  to  cover  all  the  possible  applications  of  the 
hybrid  technique  in  one  paper.  The  purpose  of  this  paper  is  to 
present  the  technique  and  demonstrate  some  of  its  facility  and  its 


accuracy.  To  accomplish  this,  various  antenna  radiators  are  considered 
in  the  presence  of  circular  cylinders.  Chapter  II  provides  the 
theoretical  background  of  the  methods  used  in  the  hybrid  technique. 

The  method  of  moments  with  the  specializations  used  is  first  described. 
Then  the  GTD  expressions  used  in  this  paper  are  presented.  In 
Chapter  III  the  hybrid  technique  is  described  in  detail.  The  very 
effective  way  in  which  GTD  is  incorporated  into  the  thin-wire  theory 
is  presented  and  demonstrated. 

Chapter  IV  is  the  results  section  in  which  the  hybrid  technique 
is  applied  to  find  the  input  impedance  of  antennas  as  a function  of 
their  distance  from  the  circular  cylinder.  For  each  case,  the  hybrid 
solution  is  checked  with  one  of  three  independent  solutions,  an 
MM-eigenfunction  solution,  image  theory,  or  experimental  measurement. 
Chapter  V concludes  the  paper  with  a summary  and  discussion. 


In  this  chapter  brief  descriptions  of  the  methods  which  are 
combined  to  form  the  hybrid  technique  used  in  this  paper  are  presented. 
Since  the  development  of  these  methods  is  not  the  subject  of  this  paper, 
the  descriptions  will  be  short  and  not  necessarily  complete.  Only 
the  parts  of  the  theories  used  in  this  paper  will  be  presented.  More 
complete  explanations  and  developments  can  be  found  in  tiie  references 
cited. 


The  moment  method  part  of  the  hybrid  technique  is  the  speciali- 
zation used  by  Richmond  in  his  thin-wire  computer  program  [3].  The 
GTD  part  of  the  hybrid  technique  is  programmed  as  presented  by 
Marhefka  [5].  These  two  references  provided  the  major  material  for 
their  respective  theories  in  the  development  of  the  hybrid  technique 
of  this  paper. 


A. 


Moment  Method 


A clear  description  of  moment  method  is  given  by  Stutzman  and 
Thiele  [6].  The  moment  method  is  a procedure  for  reducing  an  integral 
equation  of  the  first  kind  of  the  form 


I ( z 1 ) K(z.z’)  dz'  = - E1  (z) 
"over  structure 


0) 


to  a system  of  simultaneous  linear  algebraic  equations  in  terms  of 
the  unknown  current  distribution  I(z').  Once  the  current  is  known, 
determination  of  radiation  patterns  and  impedance  is  fairly  straight- 
forward. Electromagnetic  radiation  problems  can  almost  always  be 
expressed  as  integral  equations  in  the  form  of  Equation  (1).  The 
inhomogeneous  source  terms  are  on  the  right  and  the  unknown  currents 
are  within  the  integral  sign  on  the  left. 


1 . Weighted  residuals  and  the  moment  method 


A general  moment  method  procedure  can  be  accomplished  by  using 
the  method  of  weighted  residuals  [6].  The  unknown  current  I(z')  is 
expanded  on  the  structure  of  interest  using  an  appropriate  expansion 
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function,  Jn(z').  Equation  (1)  becomes  a sum  of  the  expansion  terms 
as  follows: 


N 

I 


n=l 


) K(z ,z 1 )dz' 


over  structure 


Ei(z). 


(2) 


A residual  R is  defined  to  be  the  sum  of  the  tangential  components 
of  scattered  and  incident  fields 


R = 


+ E 


tan 


(3) 


At  points  on  the  surface  of  the  perfectly  conducting  structure  of 
interest  the  residual  must  be  zero.  So 


R - l I_  Jn(z')K(z,z')dz'  + Ei(z)  = 0 


n=l 


n 

"over  structure 


(4) 


where  z is  on  the  structure's  surface.  In  the  method  of  weighted 
residuals  the  I ' s are  found  so  that  the  residual  is  forced  to  zero 
in  an  average  sense.  The  weighted  integrals  of  the  residual  are  set 
equal  to  zero  as 


Wm  Rdz  = 0,  m=l  ,2,**  *N. 


over  structure 


(5) 


where  Wr,  is  the  weighting  or  test  function.  The  choice  of  expansion 
and  weighting  functions  is  very  important  and  has  been  the  subject 
of  many  papers.  As  a rule  of  thumb,  it  is  desirable  to  choose 
expansion  functions  that  closely  resemble  the  anticipated  form  of  the 
current  on  the  structure  of  interest.  It  is  often  advantageous  to  use 
the  same  functions  for  the  weighting  functions  as  used  for  the 
expansion  functions.  This  makes  the  procedure  a Galerkin's  method. 

If  the  structure  of  interest  is  an  antenna  structure  modeled 
by  thin-wire  segments,  more  specific  expressions  can  be  described. 
Substituting  Equation  (3)  into  Equation  (5)  gives 
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E6  d*  + • E1  da  = 0, 


m = 1 ,2,*"N. 


Denoting  th^scattered  field  from  the  n-th  expansion  function  of  the 
current  by  E^,  the  mn-th  term  of  the  moment  method  impedance  matrix 
is 


W • e!  d* 

m n 


and  the  n-th  voltage  matrix  element  is 


Wm  * E1  d* 
m 


Wm  is  the  m-th  testing  function  located  interior  to  the  wire  and  on 
its  axis.  Rigorously,  the  test  function  should  be  located  on  the 
surface  of  the  wires  making  Equations  (7)  and  (8)  double  integrals  over 
the  wire  surfaces.  Placing  the  test  function  on  the  axis  is  an 
extension  of  the  electric  field  boundary  condition  for  the  sake  of 
mathematical  simplification.  This  approximation  makes  it  necessary  to 
restrict  the  method  to  wires  for  which  the  radius  is  less  than  about 
,007a. 

The  method  of  weighted  residuals  is  equivalent  to  the  steps  in 
the  usual  development  of  the  moment  method  where  a linear  operator 
and  an  appropriate  inner  product  are  defined.  In  terms  of  these  con- 
cepts, Equations  (7)  and  (8)  are 

1 = <Wm,L(J)> 

mn  m n' 

(9) 

V = <W  ,E1> 
m m’ 

where  <,>  denotes  the  inner  product  and  L is  the  linear  operator  acting 
on  Jn  to  give  E^. 


2.  Reaction  integral  equation 

A general  approach  to  the  formulation  of  boundary  value  problems 
was  presented  by  Rumsey  [7]  when  he  introduced  a physical  observable 
called  reaction.  His  formulation  resulted  in  the  reaction  integral 
equation  in  electromagnetic  theory.  Reaction  is  a measure  of  coupling 
between  two  sources.  If  the  weighting  function  is  taken  to  be  a test 
source,  then  the  impedance  ratrix  elements  given  by  Equation  (7)  may  be 
taken  as  a calculation  of  the  coupling  between  the  m-th  test  source 
and  the  field  from  the  n-th  expansion  function  or  actual  source. 
Similarly,  the  m-th  voltage  matrix  element  in  Equation  (8)  is  inter- 
preted as  the  coupling  between  the  m-th  test  source  and  the  incident 
E-field.  In  both  cases  the  m-th  test  source  current  is  reacted  with 
the  E-field  from  another  source. 

For  the  very  general  geometry  of  Figure  1,  a more  general 


Figure  1.  Test  source  inside  conducting 

body  surface  S. 


reaction  integral  equation  is  formed  from  which  Equations  (7)  and  (8) 
are  found.  Let  (Tm,Mm)  be  the  surface  current  densities  of  a test 
source  and  let  (Fm,fT,pj  be  the  fields  from  that  source.  Equivalent 
surface-current  densities  are  introduced 


Ts  = n x H" 
ct  = r x n 


(10) 


on  the  body  S.  (JS,MS)  radiate  the  fields  (E^.H5)  in  free  space  and 
replace  the  conducting  body.  The  reaction  integral  and  generalization 
of  Equation  (6)  becomes 
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J/q,^  - jVf|5>ds  * li  (Vf1  - V1'*  ■ °* 

s s 

m = 1 ,2 , " * N . (IT) 

where  the  integrals  are  over  the  test  source's  surface.  The  physical 
interpretation  is  that  there  is  zero  reaction  or  coupling  between  the 
test  source  and  the  sum  of  the  incident  and  scattered  fields.  This 
is  clearly  equivalent  to  the  weighted  residual  interpretation  of 
Equation  (5).  If  the  fields  from  the  n-th  expansion  function  of  the 
actual  source  current  are  given  by  (£p,Hp),  the  sum  of  the  N fields 
being  (F,^),  then  for  the  general  mn-th  element  in  the  impedance 
matrix 


Z = f f (J  ‘Is  - R_*H?)ds  . 
mn  J J ' m n m n' 

s 

Similarly,  the  m-th  voltage  element  is 


V = - f f (IT1  - M -TT1  )ds.  03) 

m J J m m 

s 

For  the  antenna  situations  of  interest  in  this  paper,  (F1  .TT1 ) originate 
from  impressed  currents  Xj  located  on  s. 

3.  Piecewise  sinusoidal  Galerkin  method 

The  piecewise  sinusoidal  function  as  shown  in  Figure  2 is  one  of 


Figure  2.  Typical  two  segment  dipole  with 
piecewise  sinusoidal  distribution. 
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connected  segments  and  that  each  interior  segment  contains  two  piece- 
wise  sinusoidal  functions. 

Experience  has  shown  that  the  use  of  piecewise  sinusoidal  testing 
functions  with  piecewise  sinusoidal  expansion  functions  leads  to  a 
procedure  that  is  numerically  efficient  and  highly  accurate. 

Consider  the  situation  wherein  an  arbitrary  number  of  segments  are 
coincident  with  the  z-axis.  For  example,  a dipole  antenna  with  N seg- 
ments using  Equations  (12)  and  (14)  and  the  thin-wire  assumption,  the 
mn-th  impedance  matrix  element  is 


* f l ' **  *<"■ 

• * . v ‘ <f  * . 


is  needed  to  carry  out  Equation  (15).  This  field  may  be  found  in 
a straightforward  manner  and  the  derivation  is  included  in  many 
electromagnetic  texts  including  [6].  The  geometry  is  shown  in 
Figure  4.  For  the  piecewise  sinusoidal  excitation  current  of 


z-directed  monopole  segment  source  with 
the  observation  point  (p,z). 


E = 

P 


.”,Rl  , „',R2l 


I,e  -I~e  sinhyd 


+ (I-iCoshyd-I^Je  coso^ 


I 

+ (^coshyd-^  )e  cose,. 


Ez  ■ 4»si'nM  (I,-I2cosh  d) 


+ (Ij-IjCOShyd)  ^ 


where  n is  the  impedance  of  the  medium.  The  evaluation  of  Equation 
(15)  may  be  carried  out  without  difficulty  by  numerical  integration 
or  by  Si  and  Ci  integrals  for  some  special  geometries. 


The  piecewise  sinusoidal  Galerkin  method  treated  here  is  the 
procedure  used  for  the  thin-wire  antennas  of  this  paper.  The  importance 
of  this  particular  specialization  of  the  moment  method  for  the  hybrid 
technique  will  be  discussed  in  Chapter  III. 


Geometrical  Theory  of  Diffraction 


The  geometrical  theory  of  diffraction  (GTD)  is  characterized  as  a 
high  frequency  technique  that  allows  a complicated  structure  to  be 
approximated  by  basic  shapes.  These  basic  shapes  represent  canonical 
problems  in  GTD.  Mathematically,  GTD  is  an  asymptotic  approximation 
to  an  integral  which  defines  the  electromagnetic  problem.  The  GTD 
is  a ray  optical  technique  which  allows  physical  insight  into  the 
various  scattering  and  diffraction  mechanisms  involved.  Consequently, 
the  dominant  or  significant  scattering  or  diffraction  mechanism  for  a 
given  geometry  can  be  identified  leading  to  an  accurate  engineering 
solution. 


Again,  only  the  solutions  needed  for  the  problems  solved  in 
this  paper  are  presented.  More  cases  and  derivations  including 
more  complete  explanations  are  found  in  the  literature  cited. 
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In  the  format  of  GTD,  the  total  electric  field  may  be  represented 


F = EV  + EV  + E°  . (18) 

The  F1  field  is  the  incident  field  in  free  space.  Fr  is  the  field 
reflected  from  the  surface  of  the  structure  of  interest.  F1  is  the 
diffracted  field  from  discontinuities  of  the  structure,  such  as  edges, 
corners,  or  curved  surfaces,  u1  and  ur  are  unit  step  functions  to  cut 
off  their  corresponding  field  at  the  incident  and  reflection  shadow 
boundaries.  The  extent  of  these  fields  is  determined  by  geometrical 
optics  (GO).  The  surfaces  in  this  paper  $re  all  perfectly  conducting 
and  the  medium  is  free  space.  Also  the  eJut  time  dependence  is  assumed 
and  not  explicitly  shown. 

1.  Geometrical  optics  field 

The  incident  electric  field  F1  is  considered  to  be  a spherical 
wave  for  the  equations  given  in  this  paper.  Other  cases  have  been 
treated  in  the  literature  but  the  importance  of  using  the  spherical 
wave  expressions  in  the  hybrid  technique  will  be  discussed  in  the  next 
chapter. 

The  reflected  electric  field  from  a curved  surface  S,  as  shown  in 
Figure  5,  is  given  in  geometrical  optics  terms  by 


Fr(s)  = E1(Qr)-R 


P1  p2 


(p-j+sHp^+s) 


E^Qr)  is  the  incident  field  at  the  reflection  point  (Qr)  on  the 
surface  and  R is  the  dyadic  reflection  coefficient  such  that 

= A 1 a A A 

R = eH  eH  - ^ ^ * 


Unit  vectorsAej|  and  ej|  are  parallel  to  the  plane  of  incidence  and  the 
unit  vector  ex  is  perpendicular  to  the  plane  of  incidence.  The  point 
of  reflection  Qr  is  found  from  the  laws  of  reflection  which  state 
that  the  angle  of  incidence  is  equal  to  the  angle  of  reflection.  That 


- I • n = S • n (21) 

and  that  the  incident,  reflected,  and  surface  normal  vectors  all  lie 
in  the  same  plane 


Figure  5.  Reflection  by  a smooth 
convex  surface. 


I x n = S x n . (22) 


Finding  Qr  will  be  considered  in  more  detail  in  the  last  section  of 
this  chapter. 

The  quantities  and  pj>  are  the  principal  radii  of  curvature  of 
the  reflected  wavefront  at  the  reflection  point  Qr.  Kouyoumjian  and 
Pathak  [4]  show  how  to  find  these  values  for  an  arbitrary  wavefront  by 
diagonalizing  the  curvature  matrix  for  the  reflected  wavefront  given 
by  Deschamps  [8],  The  wavefront  is  incident  on  the  curved  surface  S 
shown  in  Figure  5 at  the  reflection  point  Qr.  Unit  vectors  e-|  and 
are  in  the  principal  direction  of  S at  Qr  with  surface  radii  of 
curvature  R-|  and  R2.  For  the  case  where  the  incident  field  is 
spherical,  the  principal  radii  of  curvature  of  the  reflected  wavefront 
are  given  by  Kouyoumjian  and  Pathak  as 


S'  is  the  radius  of  curvature  of  the  incident  wavefront  at  Qr,  0]  is 
the  angle  between  the  direction  of  the  incident  ray  I and  unit  vector 
e-|,  and  02  is  the  angle  between  the  direction  of  the  incident  ray  I 
and  the  unit  vector  e2. 

2.  Diffraction  by  a curved  wedqe 

The  curved  wedge  problem  is  illustrated  in  Figures  6 and  7.  This 
report  will  consider  finite  cylinders  where  the  ends  of  the  cylinders 
are  formed  by  a flat  surface  with  a resulting  curved  wedge.  The 
diffracted  field  due  to  the  curved  wedge  is  analyzed  using  GTD 
techniques  developed  by  Kouyoumjian  and  Pathak  [4]. 

The  diffracted  field  from  the  curved  wedge  is  written  in  the 


^<s)  '^(gE)  • Be(U)  e-jks  . 


The  p is  the  distance  between  the  caustic  at  the  edge  and  the  second 
caustic  of  the  diffracted  ray.  It  is  given  by 


1 _ 1 


(I-S) 


p a sin  0 

e e 0 


where  Pg,  ne,  ae  will  be  defined  shortly.  The  diffraction  coefficients 
for  the  curved  wedge  are  extended  from  those  for  wedge  diffraction  to 
allow  the  diffracted  field  to  be  continuous  at  the  incident  and 
reflected  shadow  boundaries.  To  accomplish  this  the  appropriate 
distance  parameters  L in  each  of  the  transition  functions  make  the 
fields  continuous.  The  diffraction  coefficient  for  the  curved  wedge 
is  given  by 

-j*/4  [2  sin(^-)F[kLia(e“)] 

n /i  ii  Q \ _ e ' 


s,h  0 2n/?rcks  i ne 


0 


1 (f£) 

* cot  (tSt) 


F[kLrna+(e+)] 


F[kLroa(3+)] 


where  a(B)  = 2 cos2  6/2  and  a+(e)  = 2 cos2  (27fn-6) 


i 

i 


I t 


S'  is  the  radius  of  curvature  of  the  incident  wavefrQnt  at  Qr,  0]  is 
the  angle  between  the  direction  of  the  incident  ray  I and  unit  vector 
ei,  and  02  is  the  angle  between  the  direction  of  the  incident  ray  I 
and  the  unit  vector  e2. 

2.  Diffraction  by  a curved  wedge 

The  curved  wedge  problem  is  illustrated  in  Figures  6 and  7.  This 
report  will  consider  finite  cylinders  where  the  ends  of  the  cylinders 
are  formed  by  a flat  surface  with  a resulting  curved  wedge.  The 
diffracted  field  due  to  the  curved  wedge  is  analyzed  using  GTD 
techniques  developed  by  Kouyoumjian  and  Pathak  [4]. 

The  diffracted  field  from  the  curved  wedge  is  written  in  the 

form 

^(S)  'F'OIe)  • &E(U>  e'3ks  • <24) 

The  p is  the  distance  between  the  caustic  at  the  edge  and  the  second 
caustic  of  the  diffracted  ray.  It  is  given  by 


2 

p 


. A A , 

• (I-S) 


(25) 


where  pg,  ne,  ae  will  be  defined  shortly.  The  diffraction  coefficients 
for  the  curved  wedge  are  extended  from  those  for  wedge  diffraction  to 
allow  the  diffracted  field  to  be  continuous  at  the  incident  and 
reflected  shadow  boundaries.  To  accomplish  this  the  appropriate 
distance  parameters  L in  each  of  the  transition  functions  make  the 
fields  continuous.  The  diffraction  coefficient  for  the  curved  wedge 
is  given  by 


Dsth^**,,B0)  = 


e-3"/4 

2n/5TrRsin0 

0 


2 s1n(J)F[kLia(B“)] 
cos(J)-cos(jp) 


(26) 
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OBSERVATION 


POINT 


Figure  6.  Geometry  for  the  three-dimensional 
curved  wedge  diffraction  problem. 


Figure  7.  The  curved  wedge  in  the  plane 
perpendicular  to  e at 


I 


The  Dc  coefficient  applies  for  the  E-field  vector  parallel  to  the 
edge  (i.e. , acoustically  soft  F|Edge  = 0)  while  applies  to  the 
field  component  perpendicular  to  the  edge  (i.e.,  acoustically  hard 

=0).  The  angular  relationships  are  expressed  by  6=6 

n Edge 

The  transition  function  F(x)  is 


F(x)  = 2j  |/7|  ejx  J 

|/x| 


e"J  dz  . 


The  distance  parameter  associated  with  the  incident  field  is  given  by 


/ 1.  v i i . 2. 
i s(p  +s)p,p?sin  6 

L1  - -rA* U & . (28) 

Pg(p]+s)(p2+s) 

The  distance  parameter  associated  with  the  reflected  field  from  the 
surface  with  superscript  o (reflection  boundary  at  ) and  the 
reflected  field  from  the  surface  with  superscript  n (reflection 
boundary  (2n-l ) is  given  by 

s(p£+s)p!fp£sin2e 

Lr  - ■■r"V-L7 2 • <29> 

pgtPl+sJtPg+s) 

The  parameter  p^  is  the  radius  of  curvature  of  the  incident  wavefront 
at  the  diffraction  point  Qe  taken^in  the  plane  containing  the 
incident  ray  and  the  unit  vector  e which  is  tangent  to  the  edge  at  Q^. 

For  the  case  of  spherical  waves  p^s'.  p]  and  p^  are  the  principal 

radii  of  curvature  of  the  incident  wavefront  at  Qe-  Similarly,  p^  and 

p2  are  the  principal  radii  of  crrvature  of  the  reflected  wavefront  at 

Qe,  which  are  found  using  Equation  (23)  for  spherical  wavefront  inci- 
dence also,  the  parameter  p£  is  the  radius  of  curvature  of  the 
reflected  wavefront  at  Qe  taken  in  the  plane  containing  the  reflected 
ray  and  e.  It  is  found  using 

A A A A 

i 1 2(n*n  )( I*n) 

V - -j rt — (3°) 

pe  c€  ae  s'n  “o 
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where  n is  the  unit  vector  normal  to  the  surface  at  Qe,  ne  is  the 
associated  unit  normal  to  the  edge  directed  away  from  the  center  of 
curvature,  and  ae>0  is  the  radius  of  curvature  of  the  edge  at  Qe. 

3.  Specular  point  calculation 

The  previous  sections  presented  the  GO  and  GTD  E-field  expressions 
assuming  that  the  specular  point  Qr  and  Qe  were  already  known  on  the 
surface.  For  special  geometries  they  are  often  known  intuitively. 
Marhefka  [5]  has  presented  several  methods  for  finding  these  points 
for  more  general  geometries.  Greer  and  Burnside  [9]  have  also  investi- 
gated various  methods  of  finding  edge  diffraction  points.  The 
reflection  point  Qr  on  the  surface  of  a circular  cylinder  is  desired 
when  a nearby  source  and  field  point  are  known.  A method  presented  by 
Marhefka  was  modified  to  accomplish  this.  Figure  8 shows  the  geometry 
for  the  determination  of  the  reflection  point  on  an  elliptic  cylinder. 


( *c»  yc . *c ) 

nRSFRVATION  Pfll  NT 


Figure  8.  Geometry  for  determination  of  reflection 
point  on  elliptic  cylinder. 
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To  find  the  reflection  point  in  the  lit  zone,  the  laws  of 
reflection  as  given  by  Equations  (21)  and  (22)  are  the  starting  point. 
These  equations  are  written  such  that 

ijLlI  3 ill  = "JlI  (3i ) 

n • cT  |l|  n x 3" 


or 

(n  x T)  (rT  • (T)  + (rf  • T)  (rT  x 3)  = 0.  (32) 

Note  that  the  normalization  factors  may  be  removed.  The  incident  and 
observation  vectors  are  written  as  follows 

T(v,z)  = Ix(v)x  + Iy(v)y  + Iz(z)z 

= (a  cos  v -xs)x  + (b  sin  v -ys)y  + (z-zs)z  (33) 

and 

3(v,z)  = dx(v)x  + dy(v)y  + ^(z)  z 

= (xc-a  cos  v )x  + (yc-b  sin  v)y  + (zc-z)z  (34) 

The  normal  to  the  curved  surface  is  given  in  general  by 

/V  A A A 

n = t x z = b cos  v x + a sin  v y 
e 

or 

n = nx(v)x  + ny(v)y  . (35) 


Performing  the  necessary  dot  and  cross  products.  Equation  (32) 
becomes  two  equations 

f(v,2)  = ("xIx«yIy)(nxdy-nydx)t(nK<lxtnydy)- 


(n  I -n  I )=0 
x y y x' 


(36) 


and 
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g(v,z)  - (nxIx+nyIy)d2  + ("xdx«y<iy)Iz=0. 


(37) 


Then  the  values  of  Equations  (33),  (34),  and  (35)  are  substituted  into 
Equation  (36)  with  the  fact  a=eJv.  The  equation  is  multiplied  out 
and  like  terms  in  a are  collected  so  that  a sixth  order  polynomial  is 
obtained  given  by 


+ C^ct^  + + C i ot"^  + = 0 


(38) 


where 


C6  = (a2-b2)  [a(yc+ys)  + jb(xc+x,.)] 

c5  = - 2(a2+b2)  (x$yc+ysxc)  - j4ab[(a2-b2)  + (xsxc-ysyc)] 
C4  = a(5b2-a2)  (yc+y$)  + jb(5a2-b2)(xc+x$) 

C3  = 4(a2-b2)(xsyc+ysXc) 

C 2 = C4*  (complex  conjugate  of  C4) 

C1  = C5* 


C = 

O 6 


The  six  roots  of  the  polynomial  correspond  to  v in  that 

-1  Im(a) 

ReU) 


v = tan 


(39) 


The  value  which  is  the  true  reflection  point  is  found  by  determining 
the  v that  minimizes  the  distance  given  by  |T|  + |cT|  , which  is 
necessary  to  satisfy  Fermat's  principle.  Using  this  v parameter 
representing  the  reflection  point  in  the  x-y  plane,  the  z-coordinate 
can  be  found  from  Equation  (37)  to  be 

, ‘Vx  * nyV2s  - (V>  * Vv1  zc 

r (nxdx  + y?  ■ ‘Vx  * \V  ' 


(40) 


The  above  method  for  finding  the  reflection  point  can  be 
considered  exact.  It  is,  however,  slow  for  some  applications  and 
other  methods  are  presented  in  the  references  mentioned  above. 

In  conclusion,  the  moment  method  and  the  GTD  material  presented 
in  this  chapter  are  by  no  means  complete.  The  material  does,  however, 
provide  the  background  and  the  equations  necessary  for  the  hybrid 
technique  developed  in  the  next  chapter.  The  expressions  given  here 
are  the  ones  used  in  the  computer  programs  to  obtain  the  results  of 
Chapter  IV. 


CHAPTER  III 
THE  HYBRID  TECHNIQUE 


In  this  chapter  the  hybrid  technique  is  described  in  detail. 
The  various  problems  arising  when  GTD  is  combined  with  the  moment 
method  are  presented  and  solved.  Results  are  included  which  demon- 
strate the  effectiveness  of  these  solutions. 
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A.  Description 

The  hybrid  technique  is  used  here  to  solve  electromagnetic 
problems  in  which  antennas  or  other  discontinuities  are  located  on 
or  near  a large  conducting  structure.  The  basic  technique  was  first 
presented  in  the  literature  by  Thiele  and  Newhouse  [1],  This  method 
consists  of  casting  the  antenna  structure  in  a moment  method  format 
and  then  modifying  the  generalized  impedance  matrix  to  account  for 
the  effects  of  the  conducting  body  via  GTD. 

Following  the  notation  used  by  Thiele  and  Newhouse,  the  moment 
method  is  applied  to  the  antenna  structure  alone  by  expanding  the 
surface  current  J in  a series  of  basis  functions  J-|,  J2,  J3  *••,  such 
that 


J = 


N 

V 

L 

n=l 


*n  Jn 


(41) 


A linear  operator  L is  defined  to  relate  the  expansion  currents  to 
their  electric  fields.  A set  of  weighting  functions  W-| , W2,  W3**» 
is  selected  and  an  inner  product  is  defined  so  that. 


N 


n=l 


I_  <W  ,L(J  )>  = <Wm,E1>, 
n m n m 


(42) 


where  E is  the  field  incident  on  the  antenna.  This  is  the  m-th  row 
of  the  system  of  N equations  described  in  Section  A of  Chapter  II 
under  the  moment  method.  Equation  (42)  is  represented  as 


[Z]  (I)  = (V), 


(43) 
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The  elements  of  this  impedance  matrix  are  those  of  the  free  space 
impedance  matrix  since  only  the  antenna  structure  has  been  considered 
so  far.  These  elements  are  given  by 


mn 


• L<V>' 


(44) 


The  inner  product  forms  a unitary  space  in  which 

<J,  aE-j  + bE£>  = a<J,  + b<J,  , (45) 

where  a and  b are  complex  scalars.  If  aE]  in  Equation  (45)  represents 
L(Jn)  in  Equation  (44)  (that  is,  the  field  due  to  Jn)  and  if  bE2  in 
Equation  (45)  represents  an  additional  field  contribution  to  Zmn  (that 
is  also  due  to  Jn  but  not  due  to  the  field  arriving  directly),  then 


aE-j  + b^2> 


= <W 


m* 


L(0n) 


bL(On)> 


a=l , b=b(m,n) 


(46) 


or 


Z ' 
mn 


= <W 


m’ 


L(J  )>  + 
' n' 


<W. 


m' 


bL(Jn)> 


(47) 


The  superscript  g indicates  that  Z$n  is  added  to  each  impedance  matrix 
term  to  account  for  contributions  at  the  m-th  observation  point  due  to 
the  Jn  fields  scattered  from  the  conducting  body.  Thus  Equation  (43) 
becomes 

HZ']  [I1]  = (V)  (48) 

where  [Z 1 ] is  the  generalized  impedance  matrix  properly  modified  to 
account  for  the  presence  of  the  scattering  body  as  well  as  for  the 
antenna  itself.  Z$n  elements  are  found  with  the  aid  of  GTD.  The 
solution  of  Equation  (48)  is 

(I')  = [ZT1  (V)  (49) 


where  (I*)  is  the  current  on  the  antenna  structure  located  on  or  near 
the  conducting  body.  Thus  the  hybrid  technique  used  in  this  paper  is 
a modification  or  extension  of  the  moment  method.  This  differs  from 
other  approaches  that  consider  extension  of  GTD  via  the  method  of 
moments  [10]. 
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The  starting  point  for  the  hybrid  technique  used  in  this  paper 
is  Richmond's  thin-wire  program  [2,3].  Some  reasons  for  this  choice 
are  obvious.  It  is  a very  general  and  powerful  technique  for 
modeling  and  solving  a great  variety  of  antenna  and  scattering  problems 
with  relatively  simple  inputs  to  the  computer  code,  it  is  widely  known 
and  used,  and  its  accuracy  has  been  repeatedly  demonstrated.  The 
practical  limitation  on  its  capabilities  is  the  size  of  the  impedance 
matrix  required  by  the  geometry  of  the  particular  problem.  The 
computer  must  invert  and  solve  this  matrix.  Long  and  expensive  com- 
puter runs  become  a problem  if  the  structure  is  electrically  large. 

Another  less  obvious  reason  for  choosing  Richmond's  thin-wire 
approach  is  that  it  uses  piecewise  sinusoidal  expansion  and  weighting 
basis  functions.  The  advantages  of  these  functions  other  than  rapid 
convergence  and  the  corresponding  mode  economy,  are  made  clear  in  the 
section  of  this  chapter  on  integrating  GTD  with  the  moment  method. 


To  facilitate  a more  detailed  explanation  of  the  hybrid 
technique,  its  application  to  a specific  geometry  is  presented.  The 
case  chosen  is  a half-wave  dipole,  axially  oriented,  near  a perfectly 
conducting  circular  cylinder.  The  geometry  is  depicted  in  Figure  9. 

The  current  distribution  on  the  dipole  near  the  cylinder  is  determined 
and  from  this  the  input  impedance  is  calculated.  The  choice  of  this 
particular  geometry  and  determination  of  the  input  impedance  does  not 
imply  restrictions  on  the  hybrid  technique  described.  This  method  can 
be  applied  to  an  arbitrarily  shaped  radiator  near  any  scattering 
structure  for  which  a GTD  solution  exists.  It  may  be  used  to  find  near 
and  far  fields  and  scattering  data  as  well  as  current  distributions 
and  input  impedance. 


The  dipole  in  Figure  9 is  divided  into  segments  not  longer  than 
a quarter  wavelength  in  extent.  These  segments  are  grouped  two  at  a 
time  to  form  modes.  A moment  method  formulation  of  the  dipole  in 
free  space  is  first  carried  out  by  assuming  a piecewise  sinusoidal 
current  distribution  on  a particular  two-segment  mode.  This  test 
mode  current  generates  an  E-field  which  is  reacted  with  all  of  the 
two-segment  modes  on  the  dipole.  Each  of  these  reactions  gives  an 
impedance  matrix  term 


Rec  Mode  k 


Ik(*)  d* 


(50) 


Fj(t)  is  the  E-field  from  test  mode  j at  the  receiving  mode  k.  I|<(i.) 
is  the  expansion  current  distribution  also  assumed  to  be  piecewise 
sinusoidal  making  this  a Galerkin  method.  A row  of  the  impedance 
matrix  is  created  by  carrying  out  this  integration  for  each  expansion 
mode  on  the  dipole.  The  other  rows  are  found  by  moving  the  test 


Figure  9.  Half-wave  dipole  axially  oriented  a distance  D from  a 
perfectly  conducting  circular  cylinder  of 
radius  A. 
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current  source  to  each  mode  in  turn.  The  resulting  free  space 
impedance  matrix  is  a function  of  the  dipole  geometry  alone.  Knowing 
the  dipole  source  voltage  and  this  impedance  matrix  allows  for  the 
calculation  of  the  dipole  current  distribution  by  solving  the  set  of 
simultaneous  linear  equations.  Once  the  current  is  known,  other 
electromagnetic  quantities  follow  directly.  This  is  the  technique 
used  in  Richmond's  thin-wire  program  [3]  which  is  described  in 
ChaPter  II. 

To  calculate  the  effect  of  the  circular  cylinder,  the  hybrid 
technique  is  applied  to  find  a delta  impedance  matrix.  When  this 
delta  impedance  matrix  ([aZ])  is  added  to  the  dipole  free  space 
impedance  matrix  ([Z])  the  resulting  modified  impedance  matrix  ([Z']) 
correctly  considers  the  presence  of  the  cylinder.  Thus 

[ZMAZ]  = [Z1]  (51) 

where  [Z ' ] is  a function  of  the  dipole  and  cylinder  geometry  alone  and 
is  not  source  dependent.  Since  the  input  impedance  of  the  dipole  is 
the  quantity  to  be  found,  the  dipole  is  assumed  to  be  a radiating 
antenna  with  a known  voltage  excitation  V at  the  center  port.  This 
source  is  modeled  with  a delta  gap  voltage.  When  the  current  distri- 
bution is  determined,  the  input  impedance  then  follows  directly  from 

Z1„  ■ VZIg  <52> 

where  Ig  is  the  current  at  the  generator  port. 

The  [aZ]  matrix  is  calculated  by  the  same  method  as  [Z]  except 
that  the  E-field  which  is  reacted  with  the  expansion  modes  is  the  field 
that  is  scattered  from  the  cylinder  instead  of  the  direct  field  from 
the  test  dipole  mode.  This  scattered  field  is  the  one  resulting  from 
the  cylinder  being  illuminated  by  the  field  from  the  test  dipole 
source.  To  find  this  scattered  field  the  methods  of  GTD  are  applied. 

To  find  a particular  [AZ]  term,  the  scattered  field  must  be 
reacted  with  the  expansion  current  on  the  receiving  mode  dipole.  This 
again  is  an  integration  given  by 
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where  1 (<(*•)  is  the  piecewise  sinusoidal  distribution  on  the  receiving 
mode  dipole  k and  Fj  is  the  field  scattered  from  the  cylinder  when 
illuminated  by  the  test  source  Ij(A)  which  also  has  an  assumed  piece- 
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wise  sinusoidal  distribution.  This  integration  is  carried  out 
numerically  using  the  Newton-Cotes  equation  after  F^t)  has  been 
determined.  In  the  computer  program,  these  mode-to-mode  AZj k ' s are 
found  by  combining  the  appropriate  segment-to-segment  mutual  impedance 
terms.  Again,  this  procedure  is  modeled  after  Richmond's  thin-wire 
program.  Each  mode  impedance  term  results  from  the  superposition  of 
four  segment  impedances  making  it  possible  to  have  two  segment  "Vee" 
dipole  modes  or  dipole  modes  which  are  not  linear.  Vee  dipole  modes 
would  be  necessary  in  a more  complex  antenna  structure  than  the  dipole 
now  under  consideration.  Thus,  what  must  be  found  with  the  aid  of  GTD 
is  the  field  scattered  from  the  cylinder  at  points  along  the  receiving 
segment  resulting  from  the  test  current  on  the  source  segment. 

To  find  E5^)  using  a GTD  format  would  require  breaking  F5^) 
into  three  basic  components  a^  sketched  in  Figure  10.  F^S-)  contains 
a reflected  component  originating  from  the  source  segment  and  reflect- 
ing directly  to  the  receiving  or  observation  segment.  also 

contains  edge-diffracted  components  originating  from  .the  source 
segment  and  diffracting  off  the  cylinder  ends.  Finally 
contains  a creeping  wave  contribution  coming  from  the  source  segment, 
attaching  to  the  cylinder,  then  propagating  around  the  cylinder  as  a 
surface  wave  and  finally  shedding  to  the  observation  segment.  Other 
components  (such  as  from  the  source  to  an  edge  diffraction,  to  surface 
wave  to  another  edge  diffraction,  to  the  observation  segment)  are 
possible  but  these  contributions  would  be  minute  and  are  justifiably 
neglected. 

For  the  particular  geometry  under  consideration,  two  assumptions 
are  made.  The  cylinder  is  assumed  to  be  electrically  long  so  that 
diffraction  from  the  ends  is  negligible,  and  the  circular  cylinder  is 
assumed  to  have  a large  electrical  diameter  so  that  the  creeping  wave 
contribution  is  minute  and  may  be  ignored.  These  assumptions  are  made 
to  simplify  the  problem  used  to  describe  and  explain  the  hybrid 
technique  and  are  not  restrictions  on  solvable  geometries.  For  the 
case  in  study,  the  dominant  contribution  to  the  scattered  field  F5^) 
is  the  reflected  field  F^i.).  The  problem  has  been  reduced  to  one  of 
finding  the  reflected  field  F^)  at  observation  points  on  the 
receiving  segment  given  an  incident  field  from  the  source  segment.  To 
find  this  reflected  field,  GTD,  (or  in  this  case  GO),  is  applied. 

Referring  to  Section  B of  Chapter  II,  the  GO  field  at  an  obser- 
vation point  (as  shown  in  Figure  5)  is 
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where  {^(Qr)  is  the  field  incident  at  the  reflection  Qr  generated  by 
the  test  current  on  the  source  segment  where  R is  the  dyadic  reflection 
coefficient 
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P]  2 are  defined  in  Equation  (23).  The  field  E^Qr)  is  known  exactly 
for  a monopole  segment  with  a piecewise  sinusoidal  current  distribution. 
This  field  will  be  considered  in  detail  in  the  next  section.  The 
reflection  point  Qr  is  found  by  applying  the  technique  described  in 
Section  B of  Chapter  II.  In  applying  Equation  (54)  two  important 
points  warrant  special  attention.  First,  to  find  Qr  on  the  cylinder, 
a source  and  observation  point  must  be  specified.  The  observation 
point  presents  no  problem.  In  finding  the  segment-to-segment  impedance 
terms  a numerical  integration  is  performed  over  the  observation 
segment.  The  observation  segment  is  divided  into  a specified  number 
of  integration  sampling  points  which  also  serve  as  observation  points. 
SpecificatiQn  of  the  source  point,  however,  poses  a problem.  The 
source  of  FMQr)  is  actually  distributed  over  the  source  segment,  so 
an  assumption  must  be  made  that  the  source  appears  to  radiate  from  one 
specific  point  located  on  the  source  segment.  This  assumption  allows 
the  location  of  Qr  to  be  determined.  The  next  section  discusses  the 
actual  choice  of  the  source  point  location  and  the  ramifications  of 
this  choice. 

The  second  important  point  for  consideration  involves  the  ray 
optical  nature  of  GTD  or  GO.  The  field  ^(Or)  is  known  at  the 
reflection  point  Qr,  but  in  finding  P'(t)  only  the  components  of  FUQr) 
which  are  perpendicular  (i)  to  the  incident  ray  path  I 
are  used  as  seen  in  Equation  (54).  GTD  does  not  provide  a method  of 
including  components  of  the  field  along  the  ray  path.  The  assumption 
that  must  be  made,  then,  is  that  the  ray  path  component  of  F^Qr)  is 
negligible.  That  this  assumption  is  fairly  accurate  until  the  aipole 
gets  very  close  to  the  cylinder  is  demonstrated  by  the  results 
presented  in  the  next  section. 

Once  £*(*>)  is  known,  it  is  dotted  with  I(t)3T  on  the  observation 
or  receiving  segment.  The  segment-to-segment  impedance  terms  are 
found  by  carrying  out  the  integration  indicated  in  Equation  (53) 
numerically.  These  terms  are  combined  to  get  the  AZ  mode-to-mode 
impedance  terms,  thus  forming  the  [az]  matrix.  This  matrix,  repre- 
senting the  cylinder  effects,  is  added  to  [Z]  to  form  [Z1]  the  modified 
impedance  matrix  which  is  solved  in  the  normal  moment  method  manner. 

The  matrix  is  inverted  and  multiplied  by  the  voltage  source  column 
giving  the  desired  current  distribution.  The  input  impedance  of  the 
dipole  is  found  using  the  method  described  by  Equation  (52). 
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Combining  GTD  With  Thin-Wire  Theory 

The  hybrid  technique  as  described  in  the  previous  section  is 
applied  to  the  problem  shown  in  Figure  9 and  the  solution  of  the 
problem  is  discussed  here.  The  input  impedance  of  the  half-wave  dipole, 
axially  oriented  at  a distance  D from  a perfectly  conducting  circular 
cylinder  of  radius  A,  is  calculated.  The  real  and  imaginary  parts  of 
the  impedance  are  determined  as  a function  of  dipole  distance  D from 
the  cylinder  and  are  stored  for  plotting.  The  accuracy  of  this 
method  is  checked  by  comparing  the  plot  obtained  by  this  procedure 
with  that  obtained  by  an  independent  method. 

The  independent  method  is  a solution  by  Ersoy  and  Wang  [11]. 

The  method  consists  of  finding  the  solution  of  an  axial  dipole  near  an 
infinitely  long,  perfectly  conducting  circular  cylinder.  This  is 
accomplished  by  a technique  very  similar  to  the  hybrid  technique 
described  in  this  paper.  A delta  impedance  matrix  representing  the 
cylinder  effects  is  found  via  a moment  method  procedure.  The  method 
incorporates  the  cylindrical  Green's  function  in  the  kernel  of  the 
integral  equation.  These  eigenfunctions  account  for  the  cylinder. 

To  avoid  the  confusion  of  distinguishing  between  two  hybrid  type  tech- 
niques, this  method  will  be  called  the  "MM-eigenfunction"  technique. 

The  MM-eigenfunction  solution  is  obtained  in  the  form  of  a 
sunmation,  and  accurate  results  are  obtained  when  sufficient  terms  are 
included.  The  MM-eigenfunction  solution  is  taken  to  be  the  correct 
(exact)  one.  Figures  11a  and  b show  the  plots  for  the  case  A=A/2  and 
D is  varied  from  near  zero  to  one  wavelength.  The  agreement  between 
the  hybrid  solution  and  the  MM-eigenfunction  solution  is  good  until  the 
dipole  is  less  than  0.2A  from  the  cylinder.  The  input  impedance  of 
the  half-wave  dipole  of  diameter  0.0002A  in  free  space  with  one  piece- 
wise  sinusoidal  mode  is  73.22  + j 43.41U.  The  vertical  scales  of  the 
plots  are  in  ohms.  The  input  impedance  of  the  dipole  is  seen  to  be  a 
damped  sinusoid  oscillating  about  the  free  space  value  as  the  dipole 
moves  away  from  the  cylinder.  This  is  intuitively  logical  since  the 
energy  reflected  from  the  cylinder  to  the  dipole  should  alternately 
enhance  and  subtract  from  that  of  the  free  space  dipole,  thus  modifying 
the  dipole  current  distribution  and  giving  the  resulting  input 
impedance  plots. 

Several  other  details  of  this  problem  should  be  discussed.  The 
MM-eigenfunction  solution  used  assumed  an  infinitely  thin  dipole  whose 
current  distribution  was  exactly  cos  kz  in  free  space.  Comparing 
this  solution  with  the  hybrid  solution  required  that  the  hybrid  dipole 
also  have  exactly  the  cos  kz  distribution.  To  accomplish  this  only 
one  mode  was  used  to  model  the  A/2  dipole.  The  piecewise  sinusoidal 
current  on  this  mode  gave  the  desired  cos  kz  distribution.  This 
restriction  of  one  mode  to  model  the  dipole  was  recognized  as  a 
detraction  from  the  accuracy  of  the  hybrid  solution  and  provided 
the  motivation  for  finding  another  independent  check  as  described 


shortly.  Also,  for  these  results,  both  the  MM-eigenf unction  and 
hybrid  solutions  assumed  an  infinitely  long  cylinder.  The  hybrid 
technique  accomplished  this  by  simply  not  including  any  edge  dif- 
fraction. 

Another  important  detail,  as  pointed  out  in  the  previous  section, 
was  the  selection  of  a single  source  point  on  the  source  segment. 

Its  location  was  needed  to  find  the  field  scattered  from  the  cylinder 
in  the  segment-to-segment  impedance  term  calculations.  A logical 
choice,  and  the  one  used  here  to  find  the  curves  of  Figures  11a  and  b, 
was  the  center  of  the  source  segment. 

Three  reasons  were  postulated  for  the  breakdown  of  the  hybrid 
solution  when  the  dipole  was  closer  than  0.2A  to  the  cylinder.  The 
first,  already  mentioned,  was  that  more  modes  might  be  required  for 
better  accuracy.  That  is,  segments  shorter  than  0.25A  could  be  used  to 
model  the  dipole.  The  second  reason  was  that  the  approximation  that 
the  source  emanates  from  the  center  of  the  segment  became  less  accurate 
as  the  dipole  approached  the  cylinder.  The  third  reason  was  that  the 
ray  path  component  of  F'(Qr)  which  GO  ignores  increased  as  the  dipole 
approached  the  cylinder.  All  three  of  these  situations  would  be 
improved  if  the  dipole  was  divided  into  shorter  segments.  By  doing  this, 
the  ability  to  compare  the  results  of  the  hybrid  solution  with  those  of 
the  MM-eigenfunction  solution  is  lost  since  the  MM-eigenf unction 
solution  requires  only  one  mode. 

Another  independent  method  of  calculation  was  formulated  using 
ground  plane  image  theory.  Richmond's  thin-wires  over  an  infinite 
ground  plane  computer  code  [12]  was  the  method  used.  Consider  a 
horizontal  half  wave  dipole  at  a distance  D above  a perfectly 
conducting  infinite  ground  plane.  This  geometry  is  similar  to  that 
of  the  dipole  near  the  infinitely  long  cylinder  of  radius  A as  A gets 
very  large.  So  the  hybrid  technique  may  be  compared  with  ground  plane 
image  theory  if  a large  enough  radius  is  chosen.  In  fact,  since  the 
differences  to  be  compared  are  small,  it  is  advantageous  to  eliminate 
the  approximation  of  large  enough  A as  follows:  The  radius  of  the 
cylinder  appears  in  the  calculation  of  the  reflection  point  Qr  and  in 
the  reflected  field  F^Qr)  of  the  GO  expression.  Equation  (54).  The 
reflection  point  calculation  does  not  need  modification  for  this 
particular  geometry  since  it  will  be  independent  of  the  radius  A.  In 
the  spread  factor  part  of  Equation  (54),  and  pj?  are  the  principal 
radii  of  curvature  of  the  reflected  wavefront  at  the  reflection  point 
Qr.  They  are  given  by  Equation  (23)  for  a spherical  wavefront  incident 
and  repeated  here  as 
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The  cylinder  radius  is  R-|=A.  As  the  cylinder  radius  approaches 
infinity,  R-|-*“,  it  becomes  an  infinite  ground  plane  and  Equation  (57) 
becomes 


(58) 


Using  this  in  the  hybrid  program  is  equivalent  to  A=“.  Figure  12 
shows  the  geometry  of  this  problem. 

The  ground  plane  image  theory  solution  places  no  restriction  on 
the  number  of  modes;  modes  were  increased  to  three  and  input  impedance 
plots  were  again  determined.  As  illustrated  in  Figure  12  the  segment 
length  is  d and  the  source  point  location  on  the  segment  is  indicated 
by  x0.  For  the  particular  case  that  was  run  xo=0.5d.  Since  the  focus 
was  on  input  impedance  when  the  dipole  was  close  to  the  cylinder,  D 
varied  from  near  zero  to  only  A/2.  The  real  and  imaginary  components 
of  the  input  impedance  were  calculated  by  the  two  methods  and  the 
resulting  curves  plotted  as  Figures  13a  and  b.  The  ground  plane  image 
theory  was  considered  correct  and  is  indicated  by  the  solid  curve 
while  the  hybrid  solution  is  indicated  by  the  dotted  curve. 

These  curves  demonstrated  that  the  postulated  improvement  in 
hybrid  theory  accuracy  for  an  increased  number  of  modes  did  not 
materialize.  This  was  explained  partly  by  the  fact  that  although  the 
segments  were  smaller  and  therefore  the  errors  introduced  due  to  the 
previously  discussed  approximations  were  smaller,  they  added  up  to 
about  the  same  error  when  combined  to  include  the  entire  dipole.  To 
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Figure  12.  Geometry  of  dipole  near  cylinder  of  radius  A=°° 
showing  source  point  location  and  image 
dipole  mode. 
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Figure  13a.  Real  part  of  axial  dipole  input  impedance 
near  cylinder  of  radius  A=». 


investigate  the  effects  of  the  GTD  approximations,  the  reflected  field 
along  the  dipole  was  studied.  Specifically  the  component  tangent  to 
the  dipole,  the  one  employed  when  integrating  to  determine  the  mutual 
impedance  terms,  was  calculated  over  the  length  of  the  dipole  L.  It 
was  hoped  that  a study  of  this  aspect  of  the  problem  would  lead  to  a 
clear  understanding  of  how  GTD  combines  with  the  thin  wire  theory. 

Figure  12  also  applies  to  the  problem  of  finding  the  E-field. 
Figures  14a  and  b show  the  magnitude  and  phase  of  the  E-field  tangent 
to  the  dipole  as  a function  of  field  point  location  along  the  dipole. 

The  dipole  was  located  at  D=0.5A  from  the  cylinder.  The  center  of  the 
source  segments  was  the  location  of  the  specified  source  point  xo=0.5d. 
The  source  of  the  field  was  the  lowest  dipole  mode  on  the  halfwave 
dipole  with  segment  length  d=0.125A.  Again,  the  hybrid  solution  was 
plotted  with  dashed  lines  while  the  correct  image  theory  solution  was 
represented  by  the  solid  curve.  As  was  expected  from  the  plots  on 
Figures  13a  and  b,  for  D=.5A  the  two  solutions  agree  fairly  well. 

Figures  15a  and  b show  the  E-field  for  exactly  the  same  case  except 
that  the  dipole  was  only  D=0.125a  from  the  cylinder-,  agreement  between 
the  two  solutions  was  expected  to  be  worse  for  this  case.  From  these 
two  sets  of  figures  it  was  seen  that  the  most  severe  problem  was  in 
the  E-field  phase.  Actually,  it  was  surprising  that  the  input 
impedances  found  bv  the  hybrid  technique  were  as  good  as  shown  con- 
sidering the  E-field  errors. 

To  check  the  effect  of  the  varying  segment  sizes  used  to  model 
the  halfwave  dipole,  variations  of  the  preceding  problem  were  solved 
for  smaller  d values.  Figures  16a  and  b show  the  magnitude  and  phase 
of  the  tangential  E-field  when  d was  A/40,  or  5 times  smaller  than  the 
previous  situation.  Again,  D=0.125A,  xo=0.5d,  A=~  and  the  image 
solution  was  compared  with  the  hybrid  solution.  The  accuracy  was  not 
greatly  affected. 

In  order  to  apply  GTD  to  these  problems,  it  was  assumed  that  the 
field  at  the  observation  point  due  to  a source  monopole  segment  came 
from  a single  point  taken  to  be  the  midpoint  of  the  segment.  The  field 
was  actually  generated  from  current  distributed  over  the  source  segment. 
Since  the  current  on  the  monopole  segment  was  not  symmetric  about  the 
selected  center  point,  it  seemed  reasonable  that  a more  accurate 
source  point  could  be  chosen.  Several  alternative  source  points  were 
tried,  a logical  choice  was  the  point  of  equal  moments.  Figure  17 
shows  a piecewise  sinusoidal  distribution  sin  Bx  on  a monopole  segment 
of  length  d.  Equating  first  moments  gave  the  following  equation 

x ^ 

( m d 

j (*m-x)sin  Bx  dx  = f (x-xm)sin  Bx  dx.  (59) 
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Figure  16a.  Magnitude  of  E-field  tangent  to  axial 
dipole  along  its  length. 
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Figure  16b.  Phase  of  E-field  tangent  to  axial 
dipole  along  its  length. 
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Figure  17.  Function  for  calculation  of  equal 
moment  point  xm. 

Carrying  out  the  integration  gave  the  following  result 

x = .65966  d (60) 

m ' 


when  d=0. 125A.  From  this  new  source  point  location  the  tangential 
E-field  was  calculated  as  previously  described.  Figures  18a  and  b 
show  the  resulting  curves.  All  other  parameters  were  the  same  as 
those  for  the  solutions  presented  in  Figures  15a  and  b.  Comparing 
these  curves  shows  that  the  choice  of  source  point  location  has  a 
significant  effect  on  the  E-field  and  has  a correspondingly  sig- 
nificant effect  on  the  input  impedance.  Figures  19a  and  b demonstrate 
the  extent  of  the  effect  of  the  source  point  location.  The  dramatic 
change  in  the  curves  from  such  a slight  shift  in  source  point  location 
demonstrates  the  sensitivity  of  the  hybrid  solution  to  this  parameter. 
From  this  data  it  was  concluded  that  an  improved  method  for  integrating 
GTD  with  the  moment  method  was  essential  if  the  desired  accuracy  was 
to  be  obtained  when  the  dipole  was  located  close  to  the  cylinder. 

As  mentioned  previously,  two  important  assumptions  were  made 
in  finding  FU)  using  the  GO  expression.  To  find  Qr  a specific 
source  point  was  chosen  even  though  the  actual  source  was  distributed 
over  the  entire  segment.  E^Qr)  was  the  exact  field  at  Qr,  but  in 
determining  Er(«.)  only  the  components  of  F(Qr)  which  were  perpen- 
dicular and  parallel  to  the  incident. ray  path  were  used;  the  assumption 
was  that  the  ray  path  component  of  F(0r)  was  negligible.  A careful 
look  at  the  exact  near  field  expressions  from  a monopole  segment  with 
a piecewise  sinusoidal  current  distribution  showed  how  making  both  of 
these  assumptions  could  be  avoided. 
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Figure  18a.  Magnitude  of  E-field  tangent  to  axial 
dipole  along  its  length. 
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Figure  19b.  Imaginary  part  of  axial  dipole  input  impedance 
near  cylinder  of  radius  A=°°. 


The  fields  radiated  from  the  monopole  segments  are  given  by 
Richmond  [2].  Figure  20  shows  the  z oriented  monopole  segment.  The 
piecewise  sinusoidal  current  distribution  is 
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where  I]  or  I 2 is  zero,  y is  the  complex  propagation  constant,  d is 
the  segment  length,  and  Z]  2 are  the  endpoints.  The  resulting  fields 
at  any  point  (p,z)  are  ’ 
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where  n is  the  impedance  of  the  propagation  medium.  These  expressions 
are  for  a z oriented  segment.  They  work,  however,  for  a general 
segment  with  any  skewed  orientation  through  a simple  coordinate  trans- 
formation. To  be  correct,  these  E-fields  are  only  complete  when  added 
to  the  E-fields  of  a connecting  monopole  segment.  These  connected 
monopole  segments  make  a dipole  mode  that  is  used  in  the  thin-wire 
theory.  The  crucial  fact  to  note  about  these  E-fields  is  that  they 
may  be  separated  into  fields  emanating  from  the  two  endpoints.  More- 
over, these  separated  field  contributions  have  an  e"Y”/R  term  multiplied 
by  a pattern  factor  form,  and  are  recognized  as  spherical  waves 
emanating  from  the  endpoints.  These  observations  may  be  exploited 
with  remarkable  results  to  improve  the  integration  of  GTD  with  thin- 
wire  theory.  At  the  observation  point  on  the  dipole,  the  reflected 
field  fP^*)  will  now  be  the  superposition  of  the  contribution  from  one 
endpoint  of  the  source  segment  plus  the  contribution  from  the  other 
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endpoint.  Each  endpoint  will  have  its  own  reflection  point  on  the 
cylinder.  By  separating  the  fields  in  this  way  the  distributed  nature 
of  the  source  has  been  properly  accounted  for.  Less  obvious  is  the  fact 
that  the  ray  path  components  of  the  separated  fields  are  zero  so  that 
the  assumption  that  they  are  negligible  is  exact.  This  will  be  demon- 
strated later.  Thus  the  incident  field  from  the  source  segment  at 
the  reflection  points  is  separated  to  make  it  perfectly  compatible 
with  the  ray  optical  nature  of  GTO.  In  fact,  for  the  special 
geometry  of  the  infinite  cylinder  radius,  the  hybrid  solution  is  shown 
to  be  in  exact  agreement  with  image  theory. 

This  separation  technique  was  applied  to  the  same  case  previously 
presented  in  this  section.  The  tangent  E-field  along  the  dipole  from 
the  lowest  dipole  mode  for  the  three  mode  case  was  calculated  and 
plotted  using  the  hybrid  technique.  The  results  for  D=0.125A  are  in 
Figures  21a  and  b.  Although  it  is  difficult  to  see,  both  a solid  image 
theory  curve  and  a dashed  hybrid  solution  are  plotted  on  those  figures. 
For  this  case,  the  hybrid  solution  gives  exactly  the  correct  image 
theory  solution. 


To  analytically  justify  these  surprising  results  a careful  look 
at  the  E-field  expressions  is  required.  Working  with  the  combined 
monopole  segment  fields  of  Equations  (61)  and  (62)  is  complicated. 
Figure  22  shows  the  geometry  of  a center-fed  linear  dipole  where  the 
monopole  segments  are  considered  together.  For  a piecewise  sinusoidal 
current  distribution 
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Figure  21a.  Magnitude  of  E-field  tangent  to  axial  dipole  along 
its  length  with  endpoint  separation. 


Figure  22.  Z-directed  center  fed  dipole  source 
with  observation  point  (p,<j>,z). 

When  fields  from  two  connecting  monopole  segments  are  combined  these 
expressions  are  derived  directly  from  Equations  (61)  and  (62).  Again, 
the  fields  may  be  separated  into  contributing  fields  from  z-j,  zi 
and  z3.  Justification  of  the  separation  technique  as  well  as 
demonstration  of  zero  ray  path  contribution  can  be  shown  by  proving 
the  following  hypothesis: 


F(p,Z) 


ep  ° 


+ Ez  z = E1t1+E2t2+E3t3 


(67) 


That  is,  the  total  field  at  the  observation  point  is  postulated  to 
be  the  superposition  of  contributions  from  the  two  endpoints  and  the 
feed  point.  In  addition  the  total  field  is  said  to  be  contained  in 
components  which  are  only  transverse  to  the  ray  paths.  In  Equation 
(67)  Ep  and  Ez  are  the  known  total  fields  given  by  Equations  (65) 
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and  (66).  Let  E-j  be  the  contribution  from  point  i transverse  to  ray 
path  i (i=l,2,3).  From  the  geometry  in  Figure  22  after  some 
trigonometry 

E.  = E (R.)  sine.  - E (R • , e - ) cose.  . (68) 
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To  verify  Equation  (67),  E} , E^,  and  E3  should  give  the  total  field 
Ez  and  Ep.  Using  geometry  again 

Ez  = Ei  sinei  + E2  sin62  + E3  sin03.  (72) 

Also  using  p = Rn si ne-j  = R2Sin02  = R3Sin03  in  E],  E2»  and  E3  and  the 
trigonometric  relationship  sin^e  + cos^e  = 1,  Equation  (72)  can  be 
shown  to  be  equal  to  Equation  (66).  Similarly 


Ep  = - Eicose-j  - E2COS62  - E3COS03  (73) 

is  equal  to  Equation  (65).  Thus  Equation  (67)  is  verified  and  the 
hypothesis  is  proved  demonstrating  analytically  the  exact  agreement 
between  the  hybrid  technique  and  image  theory. 
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To  relate  this  tangential  E-field  result  back  to  the  input 
impedance  calculation,  the  case  of  a half-wave  dipole  near  the  infinite 
ground  plane  was  run  with  the  endpoint  separation  technique  incorporated 
into  the  hybrid  technique.  The  real  and  imaginary  parts  of  the  input 
impedance  were  plotted  versus  dipole  distance  from  the  cylinder  in 
Figures  23a  and  b.  The  agreement  between  the  image  and  hybrid  solution 
was  nearly  exact,  which  is  a very  significant  improvement. 

The  hybrid  technique  has  been  described  in  detail.  A method 
for  combining  GTD  with  thin-wire  theory  has  been  demonstrated  which 
integrates  the  two  so  that  the  necessary  GTD  assumptions  do  not 
hinder  the  results.  In  the  next  chapter  the  case  where  wire  antennas 
are  near  a curved  surface  will  be  examined. 


d = O.I25X 

IMAGE  THEORY 
HYBRID  TECHNIQUE 


Figure  23b.  Imaginary  part  of  axial  dipole  input  impedance 
near  cylinder  of  radius  A=». 


CHAPTER  IV 

WIRE  ANTENNAS  NEAR  CIRCULAR  CYLINDERS 
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In  this  chapter  the  hybrid  technique  will  be  used  to  solve 
electromagnetic  problems  such  as  wire  antennas  near  circular  cylinders. 
The  objective  in  this  chapter  is  to  demonstrate  the  capabilities  and 
accuracy  of  the  hybrid  technique.  To  accomplish  this  all  of  the  results 
to  be  presented  will  compare  the  hybrid  technique  solution  with  results 
obtained  for  the  same  or  similar  cases  using  an  independent  method. 

The  desire  to  compare  with  an  independent  method  is  the  main  reason 
the  particular  geometries  (e.g.,  circular  cylinder)  have  been  chosen. 

The  hybrid  technique  as  described  in  the  previous  chapter  can  be 
applied  to  much  more  complicated  geometries  than  will  be  shown  here, 
but  independent  checks  are  not  readily  available. 

Three  orthogonal  directions  or  orientations  can  be  identified  in 
relationship  to  the  cylinder.  These  orientations  are  axial  (parallel 
to  the  axis  of  the  cylinder),  radial  (perpendicular  to  the  surface  of 
the  cylinder),  and  circumferential  (tangent  to  the  surface  and 
perpendicular  to  the  axis).  Antennas  or  radiators  are  chosen  with 
orientations  to  match  one  of  these  three  directions  with  respect  to 
the  cylinder.  The  reason  for  choosing  orthogonal  orientations  to 
demonstrate  the  hybrid  technique  is  obvious.  If  the  method  can 
correctly  solve  these  three  independent  orientations  then  it  can 
solve  any  arbitrary  combination  of  them  in  the  presence  of  the 
cylinder.  Subsequently,  a general  radiator  can  be  solved  with 
confidence. 

The  electromagnetic  problems  solved  in  this  chapter  will  all 
involve  the  calculation  of  input  impedance  of  the  radiator  as  a 
function  of  distance  from  the  cylinder.  The  hybrid  method  is  capable 
of  solving  for  many  other  electromagnetic  parameters  such  as  near 
or  far  field,  current  distributions,  scattering  data,  patterns,  as 
well  as  input  impedance.  There  are  three  reasons  for  choosing  to 
find  input  impedance.  To  determine  input  impedance,  the  current 
distribution  is  first  found.  If  the  current  distribution  is  known,  all 
the  other  electromagnetic  parameters  follow  easily.  So  finding  input 
impedance  shows  that  the  other  parameters  could  also  have  been  accurately 
determined.  The  second  reason  is  that  input  impedance  is  a measurable 
quantity  and  one  which  can  be  compared  with  other  literature.  The  third 
reason  is  that  the  current  distribution  and  correspondingly  the  input 
impedance  of  the  radiators  is  very  sensitive  to  the  location  of  the 
nearby  circular  cylinder.  Thus,  solving  for  input  impedance  is  a 
good  test  of  the  capabilities  of  the  hybrid  technique. 
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Finally,  circular  cylinders  are  picked  as  the  curved  surface 
whose  presence  is  to  be  included  for  the  following  reason:  GTD  and 
therefore  the  hybrid  technique  is  capable  of  solving  many  canonical 
(and  combinations  thereof)  problems,  but,  the  circular  cylinder 
is  one  for  which  other  solutions  are  available.  The  computer  program 
developed  using  the  hybrid  technique  of  Chapter  III  and  the  theory  of 
Chapter  II  is  capable  of  solving  an  arbitrary  radiator  near  an  elliptic 
cylinder. 

In  developing  the  computer  programming  to  solve  electromagnetic 
problems  using  the  hybrid  technique,  an  effort  was  made  to  keep  the 
expressions  and  possible  geometries  as  general  as  possible.  The 
objective  was  to  develop  a program  capable  of  handling  as  broad  a 
problem  type  as  possible.  To  achieve  this,  the  thin-wire  program  of 
Richmond  [2,3]  was  chosen  as  the  starting  point  for  the  moment  method 
formulation  of  the  problem.  Some  reasons  for  that  choice  are  based 
on  the  use  of  piecewise  sinusoidal  basis  functions  as  was  discussed 
in  Chapter  II.  Another  reason  is  because  of  the  wide  familiarity, 
acceptance,  and  use  this  thin-wire  program  has.  The  hybrid  computer 
program  is  outlined  as  follows:  A wire  model  of  the  antenna  structure 
geometry  is  inputted.  The  necessary  specifications  characterizing 
the  curved  surface  are  also  read  in.  The  antenna  structure  is  sorted 
and  the  thin-wire  modes  are  set  up.  The  free  space  impedance  matrix 
is  then  found.  Next  various  appropriate  subroutines  are  called  to 
find  the  modifications  to  [Z]  due  to  the  curved  surface  starting  with 
the  GO  routine.  The  subroutines  to  find  the  modifications  use  GTD 
programs  modified  from  ones  developed  by  Burnside,  Marhefka,  Greer, 
and  others  at  the  Ohio  State  University.  The  programs  use  expressions 
from  a paper  by  Kouyoumjian  and  Pathak  [5].  The  [Z1]  matrix  is 
calculated  and  the  current  distribution  is  found.  The  hybrid  program 
was  used  on  a modified  Datacraft  computer. 

Four  specific  geometries  will  be  studied  in  this  chapter.  Three 
independent  methods  of  solution  will  be  used  to  test  the  hybrid 
solution.  These  methods  will  be  discussed  as  they  are  used.  The 
agreement  between  these  methods  and  the  hybrid  technique,  as  will  be 
shown  by  the  results  in  this  chapter,  strongly  support  the  statement 
that  an  accurate  method  for  combining  the  moment  method  treatment 
of  wire  antennas  with  the  GTD  for  curved  surfaces  has  been  found. 

A.  Axial  Dipoles 

In  this  section  the  results  of  applying  the  hybrid  technique  to 
axially  oriented  dipoles  a distance  D from  a perfectly  conducting 
circular  cylinder  of  radius  A are  presented.  The  general  geometry  of 
the  problem  is  illustrated  in  Figure  24.  The  dipole  antenna  is  half- 
wave length  in  extent  and  the  diameter  is  .0002A.  Note  that  the 
circular  cylinder  is  labeled  infinitely  long.  This  is  rigorously  true 
but  in  practice  as  long  as  the  cylinder  extends  approximately  one 
wavelength  or  more  past  the  dipole  it  appears  to  be  infinitely  long  as 
far  as  the  dipole's  current  distribution  is  concerned.  This  statement 
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Figure  24.  Half-wave  dipole  axially  oriented  a distance  D 
from  a perfectly  conducting  circular 
cylinder  of  radius  A. 


will  be  supported  by  the  results  in  the  finite  cylinder  section  of 
this  chapter. 

As  discussed  in  Chapter  III  an  MM-eigenfunction  solution  by 
Ersoy  and  Wang  [11]  is  available  for  axially  oriented  dipoles  near  an 
infinitely  long,  perfectly  conducting,  circular  cylinder.  The  solution 
is  found  by  summing  sufficient  terms  to  evaluate  an  integral  expression. 
The  integrand  contains  a Green's  function  for  the  circular  cylinder. 

The  dipole  has  a cos(kz)  current  distribution.  This  MM-eigenfunction 
solution  can  be  considered  exact  and  is  a perfect  independent  method 
to  check  the  hybrid  solution  against.  The  restriction  is  that  in  order 
to  achieve  the  cos(kz)  current  necessary  for  correct  comparison  the 
hybrid  technique  must  use  only  one  mode  or  two  segments  to  model  the 
half-wave  dipole. 

The  first  results  comparing  the  hybrid  technique  with  the  MM- 
eigenfunction  solution  are  shown  in  Figures  25a  and  b.  The  real  and 
imaginary  parts  of  the  input  impedance  of  the  dipole  were  plotted  vs. 
dipole  distance  from  the  cylinder  D.  For  this  case  D was  varied  from 
near  0 to  onehalf  wavelength  and  the  cylinder  radius  was  one  wavelength. 
The  vertical  axis  are  in  ohms,  the  solid  line  is  the  MM-eigenfunction 
solution,  and  the  dotted  curve  is  the  hybrid  solution.  The  two  methods 
give  basically  the  same  result.  The  largeness  parameter  for  the  GTD, 
or  in  this  case  the  GO,  part  of  the  hybrid  solution  is  kA.  k is  the 
wave  number  equal  to  2tt/A.  A is  the  cylinder  radius  given  in  wave- 
lengths. GTD  uses  an  asymptotic  approximation  which  is  good  when  the 
largeness  parameter  is  greater  than  one.  For  the  case  A=1.0A  shown  in 
Figures  25a  and  b the  largeness  parameter  is  about  6.3,  much  larger 
than  one  and  the  excellent  agreement  follows. 

To  demonstrate  the  effect  of  this  largeness  parameter  on  the 
hybrid  solution  three  more  cases  are  presented.  Figures  26a  and  b show 
the  input  impedance  when  A=0.5A.  Figures  27a  and  b show  the  effect 
when  A=0.25A  where  kA  is  still  greater  than  one.  Finally  Figures 
28a  and  b show  the  input  impedance  curves  calculated  by  the  two  methods 
when  kA=.785.  This  last  case  has  violated  the  largeness  parameter  con- 
straint yet  the  agreement  between  the  MM-eigenfunction  and  hybrid  solu- 
tions is  still  fairly  good.  Comparing  the  amount  of  error  in  the  hybrid 
solution  for  these  cases  as  kA  gets  smaller  it  is  seen  that  the  GO 
breaks  down  gracefully.  One  can  push  the  hybrid  technique  as  far  as  is 
consistent  with  the  desired  accuracy  of  the  particular  application  and 
not  worry  about  sudden  breakdown. 

When  the  cylinder  radius  is  very  large  the  MM-eigenfunction  solu- 
tion needs  more  terms  to  converge.  For  this  reason  and  the  need  to 
verify  the  accuracy  of  the  hybrid  solution  when  the  dipole  is  modeled 
with  more  than  one  mode,  the  method  will  be  compared  with  a ground 
plane  image  theory  solution.  In  Chapter  III  the  image  theory  method 
was  discussed.  As  A goes  to  infinity  the  cylinder  opens  to  an  infinite 
planar  conducting  surface  or  ground  plane.  Richmond  has  adopted  his 
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Figure  25b.  Imaginary  part  of  axial  dipole  input  impedance 
near  cylinder  of  radius  A=1.0x. 
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Figure  26a.  Real  part  of  axial  dipole  input  impedance 
near  cylinder  of  radius  A=0.5x. 


Figure  26b.  Imaginary  part  of  axial  dipole  input  impedance 
near  cylinder  of  radius  A=0.5A. 
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Figure  27a.  Real  part  of  axial  dipole  input  impedance 
near  cylinder  of  radius  A=0.25A. 
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Figure  28a.  Real  part  of  axial  dipole  input  impedance 
near  cylinder  of  radius  A=0.125x. 


Imaginary  part  of  axial  dipole  input  impedance 
near  cylinder  of  radius  A=0.125A. 


Figure  28b 


thin-wire  program  to  handle  this  geometry  [12],  Solving  the  case 
illustrated  in  Figure  24  when  the  cylinder  radius  was  A=10.0  gave  the 
curves  shown  in  Figures  29a  and  b.  The  solid  curve  was  calculated 
using  ground  plane  image  theory.  Again,  the  dotted  curve  is  the  hybrid 
solution.  In  each  of  these  methods  the  dipole  was  modeled  with  four 
segments  resulting  in  three  modes.  The  agreement  is  seen  to  be  nearly 
exact. 


No  matter  how  satisfying  analytical  comparisons  are,  experimental 
verification  remains  an  undisputed  and  effective  method  of  testing  a 
solution.  It  is,  as  for  the  case  of  the  square  loop  of  a later  section, 
often  the  only  independent  method  available  for  checking.  Partly  as 
further  verification  of  the  results  already  shown  in  this  section  and 
partly  as  a test  of  the  measurement  setup,  experimental  measurements 
were  made  to  determine  the  input  impedance  of  the  axial  half  wave  dipole 
located  near  a circular  cylinder.  A great  deal  of  effort  and  planning 
was  given  to  the  design  of  an  experimental  setup  with  the  hope  that 
sufficiently  accurate  measurements  could  be  achieved  in  a short  time. 

Accurate  measurement  of  the  input  impedance  of  a balanced 
feed  antenna,  such  as  a dipole,  is  very  difficult  because  the  effect 
of  the  feeding  cable  and  balancing  device  cannot  be  found  exactly  and 
thereby  be  eliminated.  To  avoid  this  problem  a monopole  over  a ground 
plane  was  used  and  equated  to  the  dipole  case  by  doubling  the  measured 
monopole  input  impedance.  The  infinitely  long  cylinder  was,  of  course, 
approximated  by  a long  cylinder  mounted  over  the  ground  plane.  Since 
impedance  as  a function  of  dipole  distance  was  desired  either  the 
monopole  or  the  cylinder  had  to  move.  The  cylinder-ground  plane 
junction  was  critically  important  and  the  monopole  was  much  smaller, 
so  the  cylinder  was  fixed  and  the  monopole  was  moved.  Knowing  tuat  the 
same  setup  was  also  to  be  used  for  the  square  loop,  the  monopole 
locations  were  carefully  chosen.  Eight  locations  were  determined  to 
be  the  maximum  number  allowable  so  that  the  adjacent  mounting  holes 
would  not  interfere  with  a measurement.  The  holes  were  located  on  a 
spiral  which  tappered  in  towards  the  cylinder  in  such  a way  as  to 
result  in  equal  clearance  angles  between  the  square  loop  locations. 

The  size  of  the  various  components  of  the  setup  as  well  as  the 
selected  frequency  of  operation  were  dictated  by  the  available  material, 
practical  mechanical  and  electrical  sizes,  and  the  available  test 
equipment.  The  four  foot  by  four  foot  ground  plane  sets  a lower  limit 
on  the  frequency  since  the  plane  had  to  be  about  two  wavelengths  or 
more  to  appear  infinite.  A five  inch  diameter  cylinder  was  the 
largest  diameter  available  and  since  a sufficiently  large  kA  was 
desired,  the  five  inches  was  picked  to  be  one-half  wavelength. 

A high  frequency  limit  of  the  experimental  setup  was  the  diameter 
of  the  monopole.  The  N-type  connector  used  to  attach  the  monopole  to 
the  ground  plane  had  a center  wire  diameter  of  one-eigth  inch.  To 
avoid  a discontinuity  right  at  the  feed  point  the  monopole  diameter 
was  chosen  also  to  be  one-eighth  inch.  Thin-wire  theory  as  used  in 
the  hybrid  technique  restricts  wire  radii  to  be  less  than  about  .007X. 
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The  frequency  of  operation  to  make  the  five  inch  cylinder  one-half  wave- 
length is  1.1811  GHz.  At  this  frequency  the  monopole  diameter  is 
.00635*.  Thus  the  dimensions  and  frequency  are  fixed.  In  order  to 
be  equivalent  to  the  half-wavelength  dipole,  the  monopole  was  cut  to 
one-quarter  wavelength  or  2.5  inches.  All  of  the  above  components  were 
carefully  machined  to  the  specified  sizes  and  assembled  to  make  good, 
smooth  electrical  contact  which  was  necessary  because  of  the  high 
operating  frequency.  The  cylinder  and  ground  plane  were  aluminum  and 
the  monopole  was  brass.  The  conductivity  of  these  materials  was 
high  enough  to  be  assumed  perfect. 

Measurement  of  input  impedance  was  accomplished  indirectly  by 
equipment  used  to  accurately  measure  the  reflection  coefficient  at  the 
monopole  port.  Figure  30  shows  a schematic  of  the  equipment  used  in 
the  measurement  setup.  Figure  31  shows  a picture  of  the  actual 
equipment  as  set  up  during  a measurement.  Note  that  the  equipment  is 
located  under  the  ground  plane  so  that  it  cannot  interfere  electrically 
with  the  measurement.  Figure  32  shows  the  various  radiators,  which 
were  used  located  about  the  cylinder. 

Again  input  impedance  was  found  by  measuring  the  complex 
reflection  coefficient  referenced  to  a short  on  the  surface  of  the 
ground  plane.  The  signal  source  was  a sweep  frequency  generator  with 
circuitry  to  lock  the  signal  level  over  a frequency  band.  The  generator 
fed  an  s-parameter  device  which  contained  a line  strecher  to  match 
the  test  line  with  a reference  line  over  a frequency  band.  The 
s-parameter  device  had  an  attached  coupler  which  compared  the  test  and 
reference  signals  and  sent  a complex  voltage  proportional  to  their 
difference  to  a signal  analyzer  with  display  units.  The  display  unit 
then  allowed  the  signal  to  be  pictured  on  a CRT  magnitude  and  phase 
plot  or  to  be  read  accurately  in  dB  on  a magnitude  and  phase  meter. 

The  measurement  procedure  went  as  follows:  the  antenna  and 
cylinder  were  secured  in  the  position  for  the  measurement.  A sweep 
frequency  signal  from  1.0  to  1.4  GHz  was  fed  to  the  antenna  through 
the  s-parameter  device.  Using  the  CRT  display,  the  line  stretcher 
was  adjusted  so  that  a short  at  the  antenna  port  appeared  as  a single 
point  at  1 .0/1 80°.  The  short  was  then  removed  and  the  reflection 
coefficient  was  observed  on  the  CRT  over  the  frequency  band  to  insure 
that  no  problems  ( i . e . , glitches  or  discontinuities)  existed  in  a 
frequency  range  about  the  frequency  of  interest,  1.1811  GHz.  The 
generator  was  then  adjusted  with  a frequency  counter  for  the  CW  signal. 
The  display  meter  replaced  the  CRT  and  the  reference  levels  were 
adjusted  while  measuring  the  short.  With  the  short  removed  the 
reflection  coefficient  was  measured  accurately  on  the  most  sensitive 
scale  settings.  Environmental  effects  such  as  unwanted  reflections 
were  checked  using  a reflecting  disk.  The  measurements  made  using 
this  procedure  and  equipment  proved  to  be  very  stable  and  repeatable. 


Figure  30.  Block  diagram  of  experimental  measurement  setup. 
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Figure  31.  Experimental  setup  showing  ground  plane  mounted  over  test  equipment. 
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Figure  32.  Various  antenna  measured  mounted  over  ground  plane  around  the  cylinder. 


To  convert  the  reflection  coefficient  R to  the  input  impedance, 
the  following  formula  was  used: 


1 - 7 1 + R 

Zin  ‘ o T 


(74) 


where  Z0  is  the  characteristic  impedance  of  the  connectors  and  cables 
attaching  the  signal  source  to  the  monopole  (Zo=50n). 


To  compare  with  experimental  measurements  several  adjustments 
to  the  dipole  data  inputted  to  the  hybrid  technique  are  necessary. 

There  are  no  modes  in  the  physical  setup  but  the  number  of  modes  used 
in  the  hybrid  technique  has  a significant  effect  on  the  input 
impedance.  To  determine  the  number  of  modes  necessary  to  get  a 
converged  number,  the  input  impedance  of  the  dipole  was  calculated  with 
D fixed  at  approximately  0.625A  and  the  number  of  modes  was  varied  from 
one  to  twenty  five.  The  resulting  plots  for  real  and  imaginary  part  of 
Zin  are  shown  in  Figures  33a  and  b.  Note  that  the  vertical  scales  have 
been  greatly  magnified  to  exaggerate  the  effect.  Choosing  9 modes  to 
model  the  dipole  will  result  in  the  real  part  being  only  0.7%  below 
the  25  mode  value  while  the  imaginary  part  will  be  within  3.0%.  The 
choice  of  9 modes  is  thus  justified. 


Another  adjustment  which  must  be  made  is  the  diameter  of  the 
dipole  to  0.0125A.  This  results  in  a length  to  radius  ratio  of  40  for 
the  dipole.  This  is  small  enough  to  be  considered  moderately  thick. 

Two  effects  may  be  significant  for  a thick  dipole.  E-field  fringing 
is  one  which  thin-wire  theory  correctly  handles  if  sufficient  modes  are 
used.  End  capacitance  is  the  second  effect  due  to  charges  on  the  ends. 
However,  one  of  the  thin-wire  assumptions  is  that  the  current  goes  to 
zero  at  the  ends  of  the  dipoles.  Forcing  the  current  to  zero  at  the 
end  incorrectly  eliminates  the  end  charge  effect.  A reasonable  approxi- 
mation to  improve  the  model  to  more  closely  fit  the  physical  conditions 
which  lets  the  current  be  non-zero  at  the  ends  is  made.  The  current  is 
allowed  to  spill  over  the  edge  of  the  end  and  go  to  zero  at  the  center. 
To  approximate  this  with  the  thin-wire  part  of  the  hybrid  program,  the 
antenna  length  is  extended  by  the  amount  of  one  radius  at  each  end. 


The  experimental  measurements  are  now  compared  with  the  hybrid 
technique  solution.  Figures  34a  and  b are  the  resulting  plots.  The 
agreement  between  the  measurement  and  the  theory  are  quite  remarkable. 


Without  discussing  why  at  this  point,  measurements  were  made  on 
a monopole  with  a radius  four  times  smaller.  This  monopole  would  be 
considered  thin  but  the  above  presented  length  adjustment  was  made 
regardless.  Its  effect  was  much  smaller  as  expected.  The  plots 
comparing  the  hybrid  technique  with  measurements  for  this  thin  dipole 
are  shown  in  Figures  35a  and  b.  Again,  the  agreement  is  very  good  with 
the  absolute  values  almost  matching  as  well  as  the  variations  matching 
quite  well. 
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Figure  35a.  Real  part  of  axial  dipole  input  impedance 
near  cylinder  of  radius  A=0.25a. 


In  this  section  then,  the  hybrid  technique  has  been  shown  to 
give  excellent  results  for  an  axially  oriented  antenna  near  a circular 
cylinder.  The  results  were  checked  with  three  independent  solutions. 
The  method  was  found  to  be  quite  accurate  with  cylinders  of  radii  from 
extremely  large  values  down  to  only  0.1 25 A. 

B.  Radial  Dipoles 


The  second  radiator  orientation  which  will  be  considered  is 
the  one  radial  to  the  cylinder.  The  antenna  will  be  a radial  dipole 
whose  center  feed  point  is  located  a distance  D from  the  infinitely 
long  perfectly  conducting  circular  cylinder.  The  geometry  of  the 
problem  is  illustrated  in  Figure  36.  The  objective  is,  again,  to 
find  the  input  impedance  of  the  halfwave  dipole  using  the  hybrid 
technique  as  0 is  varied  from  nearly  one-quarter  wavelength  to  three- 
quarters  wavelength.  In  this  case  only  one  independent  method  for 
checking  the  hybrid  solution  is  available.  This  is  the  ground  plane 
image  theory  solution  as  described  in  Chapter  II  and  also  in  the 
previous  section.  It  provides  a check  when  the  cylinder  radius  A is 
large. 

This  geometry  is  a very  special  case,  one  for  which  the  usual 
GTD  or  GO  does  not  work.  A study  of  the  geometry  as  shown  in 
Figure  36  shows  why.  Consider  any  point  on  the  dipole  as  a 
source  point  and  any  other  point,  also  on  the  dipole,  as  an  observation 
point.  The  reflection  point  on  the  cylinder  is  obvious  and  is  the 
same  independent  of  where  the  source  and  observation  points  are  on  the 
dipole.  The  reflection  point  is  where  a straight  line  extending  the 
dipole  to  the  cylinder  surface  would  touch.  The  incident  and  reflected 
ray  paths  both  lie  on  this  line  and  point  in  opposite  directions.  As 
discussed  in  some  detail  in  Chapter  II,  GTD  or  GO  does  not  consider 
field  components  along  the  ray  path  but  only  ones  transverse  to  it.  In 
fact,  fields  along  the  ray  path  are  ignored  at  the  reflection  point  by 
dotting  the  incident  field  at  the  reflection  point  with  vectors 
perpendicular  to  the  incident  ray  path.  In  finding  the  segment-to- 
segment  delta  impedance  matrix  terms,  only  the  reflected  field  com- 
ponent along  the  observation  segment  is  used.  For  this  geometry  the 
observation  segment  is  coincident  with  the  reflected  ray  path.  GTD 
or  GO  gives  only  the  components  of  the  reflected  field  transverse  to 
the  observation  segment  thus  giving  all  zero  delta  impedance  matrix 
terms.  The  overall  effect  is  that  the  hybrid  technique  using  the  GO 
as  described  thus  far  will  predict  no  change  in  input  impedance  due 
to  the  presence  of  the  circular  cylinder.  This  is,  of  course,  wrong. 

To  demonstrate  the  problem  discussed  above,  the  input  impedance 
of  the  radial  dipole  was  calculated  and  plotted  as  D varied.  The 
results  are  shown  in  Figures  37a  and  b where  the  solid  line  is  the  image 
theory  solution  and  the  dotted  line  is  the  hybrid  technique  solution. 
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Figure  37b. 


Imaginary  part  of  radial  dipole  input  impedance 
near  cylinder  of  radius  A=o°. 


As  predicted,  the  hybrid  solution  incorrectly  gives  no  effect  from 
the  cylinder  and  the  components  of  the  input  impedance  remain  constant 
at  their  free  space  values  as  the  dipole  moves. 

This  geometry  with  a radial  radiator  is  a special  case  that  can 
be  handled  by  a modification  or  extention  of  GO.  The  method  of  ex- 
tention  will  be  discussed  shortly.  The  solution  to  the  problem  then 
is  to  recognize  this  special  case  and  solve  it  separately.  The  as- 
sumption in  this  proposed  solution  is  that  the  normal  hybrid  technique 
GO  contribution  gives  an  accurate  solution  right  up  to  a delta  variation 
of  an  exactly  radially  oriented  radiator.  To  check  this  assumption  the 
case  illustrated  in  Figure  38  was  solved.  The  radial  dipole  is  tilted 
5.0°  away  from  the  exactly  radial  orientation.  0 was  again  varied  and 
the  real  and  imaginary  components  of  the  input  impedance  were  found 
using  the  hybrid  technique  and  ground  plane  image  theory.  Figures  39a 
and  b confirm  the  expected  results  that  the  hybrid  technique  holds  up 
to  a small  delta  variation  from  the  exactly  radial  orientation. 

So  that  this  special  case  can  be  handled,  the  extension  of  GO 
will  now  be  described.  Suppose  in  some  manner  it  has  been  determined 
that  an  exictly  radially  oriented  segment  is  found.  For  this  segment 
a ray  path  component  of  the  E-field  is  allowed.  How  this  field 
component  reflects  at  the  surface  of  the  cylinder  must  be  specified. 
Image  theory  shows  how  the  component  reflects  from  a planar  surface. 
Except  for  a spread  factor  due  to  the  curvature  of  the  surface,  and 
a slight  change  in  the  phase  path  length,  image  theory  is  correct.  A 
reasonable  assumption  then,  is  that  the  component  along  the  ray  path 
reflects  as  image  theory  at  the  surface  and  then  scatters  as  from  a 
plane  surface.  Since  the  ray  path  vector  lies  in  the  plane  of  in- 
cidence or  reflection  it  reflects  like  the  parallel  component  which 
also  lies  in  that  plane.  That  is  the  reflection  coefficient  is  +1.0. 
Since  the  interaction  between  a radially  directed  segment  and  the 
cylinder  is  small  the  image  theory  extension  of  GO  should  be  quite 
good  for  large  diameter  cylinders,  and  become  less  accurate  for  small 
diameter  cylinders. 

To  confirm  the  accuracy  of  this  extended  GO  methe''  the  first 
problem  considered  in  this  section  was  again  solved.  Since  all  of 
the  segments  have  exactly  a radial  orientation  so  that  the  segments 
are  col  inear  with  the  ray  path,  the  special  geometry  is  recognized. 

Ray  path  components  of  the  E-fields  are  then  included  in  the  hybrid 
program  like  parallel  components  as  described  above.  The  new 
solutions  are  shown  in  Figures  40a  and  b.  The  good  agreement  demon- 
strates the  validity  of  the  method  of  extending  GO  at  least  for  large 
A. 


A second  method  for  extending  GO  or  handling  this  special 
case  where  GO  does  not  give  answers  is  also  possible.  Again,  as- 
suming the  special  case  has  been  recognized,  the  problem  segment 
could  be  perturbed  in  orientation  to  be  slightly  off  the  ray  path 
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Figure  39a.  Real  part  of  radial  dipole  input  impedance 
near  cylinder  of  radius  A=°°. 
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direction.  The  geometry  can  now  be  handled  by  the  normal  hybrid 
technique  GO.  The  result  will  be  slightly  in  error  due  to  the 
intentional  perturbation  of  the  segment  orientation  introduced. 

This  method  has  the  advantage  that  the  curvature  of  the  cylinder 
is  treated. 

To  incorporate  this  special  handling  technique  in  the  general 
hybrid  program,  the  geometries  for  which  it  is  necessary  must  be 
identified.  The  two  cases  for  which  it  must  be  used  are  when  the 
observation  segment  is  colinear  with  a reflected  ray  path  or  when  the 
source  segment  is  colinear  with  an  incident  ray  path.  The  endpoints  of 
the  segments  and  the  reflection  points  are  known.  In  three  dimensional 
geometry  the  line  passing  through  the  two  endpoints  can  be  written  in 
parametric  form.  The  reflection  point  can  then  be  checked  with  a 
simple  logic  test  to  determine  if  it  is  colinear.  If  not,  the  normal 
GO  is  used  but  if  it  is  colinear  the  special  extended  GO  is  applied. 

Even  though  radially  oriented  antennas  present  special  problems, 
the  results  in  this  section  demonstrate  that  the  hybrid  technique  can 
approximately  solve  them.  This  is  accomplished  by  a modification  or 
extension  of  GO. 

C.  Square  Loop  Antennas 

The  last  of  the  three  orthogonal  orientations  which  must  be 
tested  and  verified  is  the  circumferential  or  phi  oriented  radiator. 

The  choice  of  antenna  with  a strong  phi  component  was  given  a great 
deal  of  thought.  It  was  desired  to  test  the  hybrid  solution  for  a 
small  cylinder  radius  as  well  as  a large  one. 

Ground  plane  image  theory  provides  the  method  of  comparison  for 
large  radii.  No  antenna  with  an  analytical  solution  was  found  when 
the  radius  was  small.  One  which  allowed  experimental  verification  was 
then  desirable.  The  square  loop  oriented  as  shown  in  Figure  41  near 
the  circular  cylinder  fits  the  bill.  The  x-y  plane  cutting  the 
geometry  gives  the  necessary  image  symmetry  to  eliminate  the  balanced 
feed  problem.  A square  half-loop  and  cylinder  can  be  mounted  on  a 
ground  plane.  The  square  loop  has  one  half  of  its  extent  almost  phi 
oriented  giving  the  required  strong  circumferential  components.  The 
square  loop  antenna  was  examined  quite  extensively  by  Richards  [14]. 
Some  of  his  results  were  used  as  a check  on  the  programs  of  this 
section  before  the  cylinder  effects  were  included. 

A look  at  the  admittance  of  the  square  loop  as  a function  of 
frequency  in  wavelengths  is  useful.  Figures  42a  and  b show  the  con- 
ductance and  susceptance  of  the  square  loop  for  a varying  frequency. 

The  wire  diameter  of  the  loop  is  one-eighth  inch  in  anticipation  of 
the  experimental  loop.  The  loop  perimeter  was  10  inches  making  it  one 
wavelength  at  1.1811  GHz.  The  vertical  scales  are  in  mi  11 -mhos. 

These  curves  agree  quite  well  with  results  by  Richards.  Maximum 
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Figure  42a.  Real  part  of  the  input  admittance  of  the 
square  loop  as  a function  of 
frequency. 


admittance  will  give  minimum  impedance  which  is  desirable  so  that 
the  antenna  will  radiate.  The  frequency  of  1.1811  GHz  is  chosen 
where  the  perimeter  is  one  wavelength.  The  loop's  dimensions  have 
now  been  specified. 

The  hybrid  technique  will  first  be  applied  to  the  geometry  of 
Figure  41  when  A is  large,  A=5.0A.  The  results  are  compared  with 
the  ground  plane  image  theory  solution  in  Figures  43a  and  b.  The 
real  and  imaginary  components  have  been  plotted  as  a function  of 
square  loop  distance  D from  the  circular  cylinder.  Note  that  the 
imaginary  impedance  is  negative  for  the  loop.  The  agreement  is  almost 
exact  verifying  the  hybrid  solution  for  large  A. 

Next,  the  results  for  A=0.25A  will  be  presented  by  comparing 
the  hybrid  solution  witn  experimental  measurement.  The  techniques  and 
procedures  described  in  section  A for  the  experimental  setup  were  used. 
The  impedance  of  the  square  loop  was  found  by  measuring  the  reflection 
coefficient  of  a square  half-loop  mounted  over  the  ground  plane. 

Figure  32  shows  a picture  of  the  loop.  The  wire  diameter  was 
one-eighth  inch  so  the  feed  point  was  continuous.  The  other  end  of 
the  loop  was  screwed  to  the  ground  plane  to  make  a good  connection. 

For  the  analytical  solution,  the  mode  stability  was  again  studied  with 
the  results  in  Table  1.  Although  it  appears  that  16  modes  are  not 
enough  for  the  real  part  to  be  accurate,  this  number  was  used  in  the 
interest  of  computer  running  time  economy. 

Table  1 


Number  of 

Modes 

Loop 

Impedance 

Modeling 

the  Loop 

Re‘zin 

) Im(Z1n 

16 

88.93 

-1  30.35 

32 

82.62 

-128.66 

48 

77.41 

-128.09 

The  results  comparing  the  hybrid  technique  with  the  experimental 
measurements  are  shown  in  Figures  44a  and  b.  The  curves  track  the 
measurements  quite  well  but  the  levels  are  shifted  making  the  number 
values  disagree.  This  level  shift  was  disturbing  when  the  remarkable 
results  of  Section  A were  recalled.  No  equivalent  reasoning  such  as 
adjusting  the  dipole  length  for  end  effects  could  be  found  for  the 
loop. 


It  was  postulated  that  the  thickness  of  the  wire  might  be  causing 
the  corners  of  the  loop  to  be  inadequately  modeled.  To  check  this, 
a loop  with  wire  thickness  four  times  smaller  was  investigated.  All 
dimensions  were  the  same  except  the  thin  loops  wire  diameter  was 
1/32".  The  admittance  curves  for  the  thin  loop  case  are  shown  in 
Figures  45a  and  b.  Note  that  1.1811  GHz  is  still  a good  frequency 
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Real  part  of  square  loop  input  impedance 
near  cylinder  of  radius  A=5.0A. 
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Imaginary  part  of  square  loop  input  impedance 
near  cylinder  of  radius  A=5.0x. 


Figure  43b 


1200 

FREQ  ( MHz  ) 


1700 


Figure  45a.  Real  part  of  the  input  admittance  of  the 
square  loop  as  a function  of  frequency. 


choice  and  resonance  is  much  sharper  than  for  the  thick  loop.  Figure 
46  shows  a picture  of  the  thin  loop.  The  end  opposite  the  feed  point 
is  soldered  to  a brass  screw  set  in  the  ground  plane.  Mode  stability 
was  again  studied  with  the  results  in  Table  2.  The  sixteen  modes  for 
the  thin  loop  is  more  accurate  than  for  the  thick  loop. 


Table  2 


Number  of  Modes 

Input  Impedance 

Modeling  the  Loop 

Re(zin) 

lm(Z.  ) 
v m 

16 

101.255 

-139.802 

32 

99.379 

-139.315 

48 

98.454 

-139.199 

The  results  for  the  thin  loop  case,  comparing  the  hybrid 
technique  with  experimental  measurements  are  shown  in  Figures  47a  and 
b.  The  postulated  improvement  did  not  materialize.  The  curves,  again, 
track  the  points  well  but  a level  shift  is  still  apparent. 

During  the  measurement  process  it  was  observed  that  a slight 
change  in  frequency  resulted  in  a significant  change  in  the  input 
impedance.  This  follows  logically  from  the  admittance  curves  where 
it  can  be  seen  that  the  operating  frequency  is  at  a rapidly  varying 
point.  The  inability  to  get  exact  level  agreement  is  at  least  partly 
explained  by  the  sensitivity  at  the  frequency  of  operation.  An 
attempt  was  made  to  avoid  this  problem  by  carrying  out  the  measurements 
at  a frequency  low  enough  to  be  away  from  the  sensitivities  of 
resonance.  A frequency  of  885.8  MHz  was  chosen  which  can  be  seen  from 
Figures45a  and  b to  be  well  off  resonance.  The  results  are  shown  in 
Figures  48a  and  b.  Agreement  is  quite  erratic.  The  problem  is  that 
Z-jn  for  this  off  resonance  case  is  very  large  compared  with  the  Zo=50ft 
of  the  connecting  cable.  The  mismatch  causes  much  of  the  energy  to 
reflect  back  down  the  line  and  very  little  radiates.  Small  radiation 
means  weak  interaction  with  the  cylinder.  The  measurement  environment 
is  too  noisy  and  the  equipment  is  too  inaccurate  to  measure  this  weak 
interaction  properly. 

In  an  attempt  to  further  investigate  the  problems  around 
resonance  by  measuring  free  space  loop's  input  impedance  over  a fre- 
quency band  the  real  problem  was  discovered.  The  section  of  line 
connecting  the  antenna  to  the  s-parameter  device  is  assumed  to  be  ideal. 
Only  its  length  is  compensated  for  by  setting  the  reference  to  a 
short  at  the  antenna  port.  It  was  found,  however,  that  this  section 
of  line  which  was  constructed  of  a group  of  connectors  was  not  ideal. 
Different  connectors  of  the  same  type  gave  different  amounts  of  phase 
shift.  The  amount  of  change  in  phase  was  also  dependent  on  how  much 
the  antenna  was  radiating.  Correspondingly,  it  was  frequency  de- 
pendent. Physically,  the  difference  probably  resulted  from  the  dis- 
continuities where  the  connectors  touched  being  unique  to  each  connector. 
The  solution  of  the  problem  would  involve  using  a more  ideal  connector 
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Figure  47a.  Real  part  of  square  loop  input  impedance 
near  cylinder  of  radius  A-0.25A. 
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Imaginary  part  of  square  loop  input  impedance 
near  cylinder  of  radius  A=0.25x. 
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Figure  48b.  Imaginary  part  of  square  loop  input  impedance 
near  cylinder  of  radius  A=0.25A. 


in  this  critical  part  of  the  setup.  Remeasurement  was  not  necessary 
here  because  the  results  already  show  that  the  hybrid  technique  works 
very  well  and  has  been  experimentally  verified. 

Thus,  phi  oriented  radiators  have  been  checked  for  both  large  and 
small  cylinder  radii.  The  hybrid  technique  has  now  been  demonstrated 
to  be  accurate  for  all  three  orthogonal  orientations.  Its  use  on 
arbitrary  radiators  near  the  circular  cylinder  can  now  proceed  with 
confidence. 

D.  Finite  Length  Circular  Cylinders 

In  this  section  finite  length  cylinders  will  be  handled  using 
the  hybrid  technique.  To  account  for  the  finite  cylinder  length,  the 
GTD  for  curved  edges  is  used  as  described  in  Chapter  II  to  find 
delta  impedance  matrices  as  described  in  Chapter  III  which  properly 
modify  the  free  space  matrix  to  include  the  cylinder  end  effects.  The 
formulation  of  the  part  of  the  hybrid  technique  which  finds  the  cylinder 
end  effects  is  modeled  after  the  part  which  finds  the  reflection  effect. 
The  formulation  was  described  in  detail  in  Chapter  III. 

To  test  the  finite  cylinder  capability  of  the  hybrid  technique, 
it  was  applied  to  find  the  input  impedance  of  an  axially  oriented 
dipole.  One  reason  for  picking  this  radiator  orientation  was  that 
the  diffraction  points  on  the  cylinder's  edges  were  known  directly. 

The  problem's  geometry  is  illustrated  in  Figure  49.  The  diffraction 
points  are  on  the  cylinder's  end  edges  in  the  x-z  plane  with  x positive. 
The  cylinder  half  height  is  CH,  its  radius  is  A,  and  the  half-wave 
dipole  is  a distance  D from  the  cylinder.  For  the  cases  in  this 
section,  A will  be  fixed  at  one-quarter  wavelength. 

The  only  analytical  method  available  for  solving  a finite 
cylinder,  other  than  the  hybrid  technique,  is  the  moment  method. 

Computer  limits  on  storage  and  running  time  prohibit  its  use  on  a 
cylinder  as  electrically  large  as  the  case  to  be  solved.  So  no  inde- 
pendent method  is  available  to  check  the  hybrid  solution  for  this 
geometry.  Experimental  measurement  will  prove  impossible  because 
the  effects  of  the  finite  ends  will  be  too  small  to  measure  accurately. 
The  method  chosen  to  verify  the  results  was  to  compare  the  hybrid 
method  solution  of  the  finite  length  cylinder  case  with  the  hybrid 
solution  of  the  infinite  length  cylinder.  Intuitive  observations  will 
then  be  made. 

The  first  results  to  be  presented  are  for  a cylinder  height 
CH=0.375x  where  D is  varied  from  near  zero  to  one  wavelength.  The 
input  impedance  plots  are  given  in  Figures  50a  and  b.  The  solid  curves 
are  for  an  infinite  length  cylinder.  Chopping  the  cylinder  off 
creates  a small  perturbation  in  the  input  impedance  of  the  dipole. 

This  result  seems  reasonable. 
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Figure  50a.  Real  part  of  axial  dipole  input  impedance 
near  cylinder  of  radius  A=0.25A. 


Another  approach  was  tried  to  help  verify  the  results.  The 
dipole  distance  from  the  cylinder  was  held  fixed  and  the  cylinder 
height  CH  was  varied  from  just  above  the  dipole  to  one  wavelength 
higher.  Figures  51a  and  b show  the  results  for  D=0.3A.  The  solid 
line  is  the  finite  cylinder  case  here.  Figures  52a  and  b show  the 
case  D=0.6A  results.  These  two  sets  of  figures  are  particularly 
convincing  in  that  the  effects  of  the  cylinder  ends  fade  to  zero  as 
the  cylinder  gets  longer  as  it  must. 

Finally,  from  those  curves  another  cylinder  height  CH=.45A  was 
picked  and  the  input  impedance  was  again  found  as  D was  varied.  The 
results  are  plotted  in  Figures  53a  and  b.  Although  the  hybrid  finite 
cylinder  solution  has  not  been  verified  by  an  independent  method, 
the  reasonable  results  of  this  section  lend  confidence  to  the  solution. 


CHAPTER  V 
DISCUSSION 


A hybrid  technique  has  been  presented  which  solves  electromag- 
netic problems  in  which  an  antenna  or  scatterer  is  located  on  or  near 
a conducting  body.  The  general  technique  has  been  applied  to  find 
the  input  impedance  of  various  radiators  located  near  a perfectly 
conducting  circular  cylinder.  The  application  of  the  hybrid  technique 
to  these  specific  problems  in  no  way  implies  restrictions  on  its  use. 

The  method  can  be  equally  well  used  to  find  near  and  far  fields, 
current  distributions,  and  scattering  data  for  problems  involving 
conducting  bodies  of  rather  arbitrary  shape. 

The  hybrid  technique  in  this  paper  used  the  format  of  a powerful 
thin-wire  computer  program  to  model  the  antenna  structure.  The  free 
space  impedance  matrix  was  found  using  the  normal  moment  method 
technique.  Then  the  effects  of  the  conducting  body  were  incorporated 
into  delta  impedance  matricies.  A delta  impedance  matrix  was  found 
for  each  type  of  contributing  scatter  from  the  conducting  body,  such 
as  reflection,  edge  diffraction,  etc.  The  delta  impedance  matricies 
were  added  to  the  free  space  matrix  to  account  for  the  conducting 
body.  The  modified  impedance  matrix  was  then  used  in  the  normal  moment 
method  way  to  find  the  current  distribution  and  input  impedance  of  the 
antenna. 

The  delta  impedance  matrix  terms  were  found  by  reacting  an 
expansion  current  with  the  E-field  scattered  from  the  conducting  body 
due  to  a test  current  source.  Piecewise  sinusoidal  expansion  and  testing 
functions  were  used.  The  GTD  E-field  at  an  observation  or  integration 
point  on  the  receiving  mode  was  composed  of  the  superposition  of  end 
and  feed  point  contributions  from  the  source  mode.  Separating  the 
contributions  in  this  way  accounted  for  the  distributed  nature  of  the 
source  and  resulted  in  spherical  waves  emanating  from  the  points  which 
were  completely  compatible  with  the  ray  optical  nature  of  GTD.  This 
way  of  combining  thin-wire  theory  with  GTD  gave  results  that  were  far 
more  accurate  than  other  formats.  Analytical  justification  as  well 
as  verification  by  example  were  presented  in  Chapter  III. 

The  accuracy  and  facility  of  the  hybrid  technique  were  shown  in 
Chapter  IV  by  solving  for  the  input  impedance  of  various  radiators  as 
a function  of  their  distance  from  a circular  cylinder.  Three 
orthogonal  orientations  were  identified  and  antennas  to  match  them  were 
used.  With  the  hybrid  technique  verified  for  these  orthogonal  orien- 
tations, clearly  a more  general  radiator  with  an  arbitrary  combination 
of  orientations  could  be  solved  with  confidence. 
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An  axially  oriented  half-wave  dipole  near  the  circular  cylinder 
was  solved  by  the  hybrid  technique.  Wlien  the  cylinder  radius  was 
large  the  hybrid  technique  was  compared  with  an  image  theory  solution. 

For  small  cylinder  radii  the  hybrid  solution  was  compared  to  an  exact 
MM-eigenfunction  solution.  In  both  cases  the  agreement  was  excellent. 

For  this  orientation  the  solution  was  also  compared  with  careful  ex- 
perimental measurements.  The  measurements  were  accomplished  by 
doubling  the  input  impedance  of  a quarter-wave  monopole  mounted  with 
the  cylinder  over  a ground  and  measuring  the  reflection  coefficient. 
Results  for  both  a thick  and  a thin  dipole  agreed  very  well  with  the 
measurements. 

A radially  oriented  dipole  presented  a special  problem  in  that 
ordinary  GO  predicted  no  effect  on  input  impedance  when  the  dipole  was 
close  to  the  cylinder.  This  incorrect  solution  was  due  to  the  non-zero 
ray  path  E-field  which  GO  ignores.  GO  was  extended  to  include  a ray 
path  field  for  this  special  orientation.  The  component  reflects  like 
image  theory  at  the  surface  and  then  scatters  as  from  a plane  surface. 
Tilting  the  problem  segment  slightly  from  its  exactly  radial  orientation 
offered  another  solution.  Comparing  the  hybrid  solution  extended  in 
these  ways  with  image  theory  verified  the  techniques. 

The  circumferential  orientation  was  verified  with  a square  loop 
lying  tangential  to  the  circular  cylinder.  For  large  cylinder  radii 
the  hybrid  solution  agreed  almost  exactly  with  image  theory.  For 
smaller  cylinder  radii  the  solution  was  compared  with  experimental 
measurements  made  on  a half-loop  mounted  with  the  cylinder  over  the 
ground  plane.  Both  the  thick  and  thin  loop  solutions  tracked  the 
impedance  variations  measured.  The  absolute  levels  were  shifted 
slightly  for  two  reasons.  The  frequency  of  operation  was  near  resonance 
where  rapidly  varying  impedance  makes  level  matching  difficult.  Also, 
a section  of  line  which  was  assumed  to  be  ideal  had  a frequency  dependent 
phase  shift. 

Thus,  three  orthogonal  orientations  were  verified  and  the  hybrid 
technique  was  clearly  shown  to  be  very  accurate.  A finite  length 
cylinder  was  also  investigated  with  the  hybrid  technique  when  an  axial 
dipole  was  nearby.  The  results  were  compared  with  the  infinitely  long 
cylinder  case.  As  expected,  the  finite  length  effects  died  out  as 
the  cylinder  became  longer. 

Several  minor  limitations  and  restrictions  on  the  hybrid  technique 
can  be  pointed  out.  The  antenna  structure  near  the  conducting  body 
is  limited  in  electrical  size  since  the  antenna  must  be  modeled  by 
wire  segments  no  more  than  one-quarter  wavelength  in  extent.  The 
number  of  segments  determines  the  number  of  modes  and  correspondingly 
the  size  of  the  impedance  matrix  the  computer  must  work  with.  The 
conducting  body  must  be  a canonical  shape  or  combination  thereof  for 
which  GTD  solutions  exist  for  the  canonical  parts.  Besides  the  field 
expressions,  a method  of  finding  the  specular  points  on  the  surface 
must  also  be  available. 
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In  spite  of  the  restrictions  discussed,  a very  powerful  hybrid 
technique  has  been  presented  for  combining  the  moment  method  treatment 
of  wire  antennas  with  the  GTD  for  curved  surfaces.  The  method  has 
been  demonstrated  to  be  accurate  and  versatile. 
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